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Shock  Wave  Generation  in  Water 
for  Biological  Studies 

Dongkoo  Kang,  Jong  Bong  Nah,  Michael  Cho,  and  Shu  Xiao,  Senior  Member,  IEEE 


Abstract — We  have  used  a  pulsed  power  generator  based  on 
stacking  MOSFET  switches  to  create  breakdown  in  water  and 
generate  shock  waves.  The  pulses  with  an  exponential  decay 
have  a  maximum  voltage  of  4  kV.  With  the  pulse  application 
of  approximately  50  the  breakdown  voltage  is  approximately 
3  kV,  Postbreakdown  events,  such  as  shock  wave  emission  and 
bubble  expansion  and  collapse  were  recorded  by  a  fast  camera. 
The  shock  speed  was  recorded  as  1.5  km/s  and  the  maximum 
bubble  diameter  is  2.3  mm.  This  shock  wave  generator  can  be 
used  on  the  microscope  and  provides  a  means  of  studying  the 
shock  impact  on  biological  cells. 

Index  Terms — Biological  studies,  shock  wave,  stacking 
MOSFETs,  water  breakdown. 

I.  Introduction 

HOCK  waves  have  been  used  in  dismantling  kidney  stones 
(lithotripsy)  [1],  tumor  destruction  [2J,  and  drug  or  gene 
delivery  [3],  [4).  The  high  gradient  of  the  pressure  carried 
by  shock  waves  causes  complicate  reflections  and  refrac¬ 
tions  when  they  impinge  sonic  impedance  discontinuities. 
A  rarefaction  wave  follows  the  compression  wave  and  carries 
a  negative  pressure,  causing  the  coexistence  of  compression 
and  expansion.  It  is  such  compressive  and  tensile  forces 
that  cause  most  of  the  damages.  In  the  liquid  environment, 
such  transition  of  positive  to  negative  pressure  could  induce 
microbubbles  and  microstream  force,  which  are  responsible 
for  the  pore  formation  of  cells  (sonoporation)  [51, 

Generating  a  high  gradient  pressure  front  can  be  done 
by  shock  tubes  or  gas  guns.  They  are  also  used  in  gen¬ 
erating  a  long  lasting  blasting  wave  up  to  milliseconds  or 
seconds.  In  the  liquid  environment,  for  example,  under  water, 
a  piezoelectric  transducer,  or  an  arc  discharge  can  be  used  for 
repeatable  use  in  a  controllable  condition  [6].  For  the  purpose 
of  studying  the  impact  of  pressure  waves  to  biological  cells, 
we  have  used  the  arc  discharge  considering  the  fact  that  it  is 
possible  to  construct  a  low-profile  shock  wave  system  on  the 
microscope  stage.  This  would  be  of  an  interest  in  studying 
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the  response  of  biological  cells  to  shock  waves,  in  particular, 
cells  cultured  on  a  coverslip.  Immediate  effects  can  be  studied 
by  performing  real-time  measurement  on  the  microscope  and 
therefore  such  system  will  be  useful  in  understanding  the 
interaction  of  shock  waves  with  biological  cells. 

Shock  waves  generated  by  electric  discharges  in  the  liquid 
were  studied  as  a  byproduct  of  high  power  switching  using 
liquid  as  medium  [7]-[9],  Numerous  data  have  been  generated 
and  useful  conclusions  have  been  drawn.  Most  of  the  data 
were  focused  on  high  voltages  {above  10  kV)  and  high  peak 
currents  (up  to  kiloampere)  [IO]-[!2].  We  instead  use  lower 
voltages  (4  kV)  and  low  currents  {<50  A)  to  produce  shock 
waves  in  a  culture  dish.  The  pulsed  power  source  was  based 
on  a  capacitor  discharge  switched  by  four  MOSFET  switches 
stacking  in  series.  Our  aim  was  to  establish  a  consistent, 
reliable  shock  wave  source.  Specifically,  we  have  studied  the 
temporal  evolution  of  shock  waves  and  bubble  expansion 
through  a  high  speed  camera.  The  pressure  change  over 
distance,  the  bubble  radius  over  distance,  and  their  relation 
at  the  onset  of  the  shock  wave  are  examined.  We  estimate  the 
shock  pressure  from  the  speed  of  the  shock  wave  and  bubble 
expansion. 

IL  Theoretical  Considerations 

Using  pulsed  power  systems  to  generate  shock  waves 
in  water  relies  on  electric  breakdown.  One  may  also  pro¬ 
duce  pressure  waves  when  streamers  are  initiated  1 1 3 J. 
Higher  pressure  is,  however,  produced  when  electric  break¬ 
down  occurs  and  the  electric  field  can  be  as  high  as 
I  MV/cm  [141,  [15].  Rapid  heating  occurs  when  the  current 
increases  and  the  voltage  across  the  electrode  gap  collapses. 
Significant  amount  of  electric  energy  deposits  in  the  plasma 
channel  in  the  form  of  heating  and  vaporizing  the  liquid 
by  increasing  the  internal  and  kinetic  energies.  The  rapid 
expansion  of  the  plasma  channel  passes  momentum  to  the 
surrounding  water,  compressing  the  water,  and  forming  shock 
waves.  During  earlier  time  of  the  expansion  of  the  plasma 
channel,  the  pressure  behind  the  shock  wave  front  is  smaller 
than  the  pressure  at  the  boundary  of  the  plasma  channel  .  As  the 
shock  waves  detach  from  the  plasma  channel,  they  propagate  at 
a  supersonic  speed  very  near  the  electrodes  but  become  sonic 
as  the  distance  increases.  Even  at  the  sound  speed,  the  wave 
could  still  carry  a  pressure  greater  than  the  ambient  pressure. 
While  the  shock  wave  emission  occurs  much  faster,  the  plasma 
channel  or  bubble  left  behind  still  has  significant  pressure  and 
keeps  growing  in  volume  until  it  reaches  its  peak  size  before 
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collapsing.  A  second  shock  wave  may  be  observed  after  the 
collapse  of  the  gas  bubble  [10]-[16], 


A.  Bubble  Dynamics  and  Pressure 


In  the  stage  of  bubble  expansion,  the  volume  grows  from 
streamer  channels,  which  is  created  in  the  prebreakdown 
phase.  The  energy  balance  equation  is  1 17] 


p~\ 


PcV<  + 


1 


Vcpc  =  i2r, 


(1) 


where  P<  is  the  pressure  of  the  plasma  channel,  /i  is  the 
effective  adiabatic  exponent,  and  Rv  is  the  measurable  resis¬ 
tance  of  the  plasma  channel  and  depends  on  time.  A  more 
complicated  form  of  the  plasma  channel  resistance  including 
nonunitorm  density  and  temperature  distribution  can  be  found 
in  [ IS]*  Here,  we  will  not  consider  such  case.  Rather,  we 
obtain  its  expression  from  measurement.  Vc  is  the  volume  of 
the  spherical  plasma  Vt  =  4jrr|/3.  The  largest  heating  occurs 
in  the  resistive  phase  of  the  plasma  channel. 

The  current  allows  the  connection  of  the  plasma  channel 
resistance  with  the  external  circuit  that  provides  the  energy 

LI  +  Ri  +  i  /  t  dz^Uo  (2) 

L  J  n 

where  R  is  the  total  resistance  of  the  circuit,  R  =  Rk  +  Rc. 
Rk  is  the  resistance  placed  in  series  with  the  gap  and  in  our 
case  is  100  Q.  The  plasma  channel  is  called  bubble  in  its  later 
phase.  The  bubble  motion  started  from  rest  is  irrotational  and 
the  equation  of  continuity  gives 


V  *  w  ^  0  (3) 

where  ;7  is  the  velocity  vector  of  the  bubble.  Once  a  velocity 
potential,  is  introduced:  u  =  V$,  the  continuity  equation 
becomes  Laplace’s  equation 

V24>  =  0.  (4) 


In  general,  with  Green's  functions,  we  can  solve  this  veloc¬ 
ity  potential  at  any  point  (r)  inside  the  bubble  using  Green's 
theorem,  which  allows  the  differential  equation  to  be  converted 
to  an  integral  equation 

(?)  =  ( 1  /4k)  fQ(Gd<p/dn  -  <f>dG/dn)dS.  (5) 

The  integral  is  performed  over  the  boundary  £1,  provided  the 
boundary  conditions  are  known.  They  typically  include  rigid 
boundaries  (no  flow  condition:  V0  ■  n  —  0,  n  is  the  boundary 
normal)  and  a  free  boundary  on  the  bubble  surface.  Using 
the  momentum  equation,  the  free  boundary  condition  can  be 
written 

Btp/dt  +  1  /2|V0|2  +  ip,.  +  pr)/p  +  gZ=  Pca/P-  (6) 


Fig.  I .  illustration  of  Green's  functions  for  a  free  field  source  and  a  source  on 
top  of  a  rigid  boundary,  (a)  In  a  spherical  expansion,  the  sources  on  the  bubble 
surface  form  a  symmetry  from  the  upper  half-plane  to  the  lower  half-plane. 
<  b>  Bubble  expansion  on  the  surface  of  a  boundary,  (c)  Effect  of  boundary  can 
be  considered  hy  adding  symmetrical  images  in  the  lower  half  plane,  which 
resembles  (,a|. 


coverslip.  The  bubble  dynamics  share  the  same  characteristics 
as  that  for  the  discharges  in  free  held.  This  is  so  because  their 
Green’s  functions  are  the  same. 

For  the  bubble  in  the  free  field,  the  Green's  function  for  the 
source  on  the  bubble  surface  is 

G  =  l/\r-r\  +  ]/\r-r"\  (7) 

where  r  locates  in  the  upper  half-plane  and  r  in  the  lower 
half-plane.  For  the  bubble  generated  on  the  coverslip  surface, 
the  coverslip  can  be  considered  by  the  method  of  images  [16], 
The  sources  like  r  ,  and  their  images,  like  r  with  an  equal 
strength  create  the  same  boundary  condition  as  the  rigid 
boundary  {Fig.  I).  The  corresponding  Green's  function  is 
identical  to  (7).  Therefore,  we  can  drop  the  integral  over  the 
rigid  boundary,  but  integrate  over  a  full  sphere  surface  in  the 
free  field.  Therefore,  the  solutions  for  both  cases  should  be 
the  same  based  on  (5)  and  (7),  Such  conclusion  is  also  drawn 
in  [20].  A  more  concise  expression  for  a  bubble  growing  in 
the  free  field  (also  applied  in  our  case)  is  developed  according 
to  the  square  root  of  time  [21] 

rc(t)  =  an.  (8) 

In  this  simplified  bubble  expansion  relation,  we  can  estimate 
the  coefficient  a  from  the  images  of  bubble  expansion.  We  note 
that  in  (8),  as  t  —  0,  the  radius  of  the  bubble  becomes 
zero,  which  implies  that  the  transition  of  streamer  to  arc 
is  instantaneous.  One  can  also  solve  a  more  complicated 
form  of  bubble  radius  [11],  [22],  which  is  a  function  of  the 
enthalpy  difference  between  the  liquid  at  the  bubble  wall  and 
at  hydrostatic  pressure. 

Assuming  the  decrease  in  resistance  of  the  arc  is  due  to  the 
increase  of  conductivity  and  the  expansion  of  the  volume  of 
the  arc,  (8)  leads  to  the  resistance  of  the  plasma  channel  [23] 

RAt)~bftl.  (9) 


Equation  (6)  is  also  called  Bernoulli  equation,  which  allows 
<j>  to  be  determined  at  later  time  t  +  <5f,  so  the  bubble  solution 
through  time  can  be  obtained  through  iteration  [19], 

In  this  paper,  we  will  not  use  this  approach  to  solve  the 
bubble  dynamics,  but  the  outlined  process  above  points  to  an 
important  conclusion  relevant  to  our  study.  That  is,  in  our 
experiment,  we  created  discharges  on  the  surface  of  a  glass 


The  constant  b  can  be  estimated  from  the  resistance  of 
the  water  arc  by  measuring  the  voltage  and  current  across 
the  electrodes.  Equations  (1),  (2),  (8),  and  (9)  determine  the 
pressure  development  as  the  energy  from  the  circuit  deposits 
in  the  form  of  Joule  heating. 

The  initial  conditions  for  the  unknowns  are  given:  /  =  0, 
dljdt  = Uq/L ,  Pt.  ^  100  kPa,  and  Vv  =  0, 
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In  solving  the  equations,  the  parameters  need  to  be 
normalized 

l=7rAJ=7W^RMZ’f  =  F/nc° 

where  po  is  the  unperturbed  water  density  and  cq  is  the  sonic 
velocity  in  water. 


B.  Shock  Wave  Pressure 

The  pressure  and  particle  velocity  at  the  gas- water  boundary 
are  equal  at  all  limes  1 16].  Knowing  the  pressure  of  the 
expansion  bubble,  its  radius  and  speed,  one  can  obtain  the 
shock  pressure  at  a  distance  /■  [  1 1 J 

-B  (10) 


p(r)  =  A 


n  +  l  n  +  I 


V  rcn  / 


where  G  =  /?/,(//  +  U2/ 2}  is  a  constant  along  a  path  or 
a  characteristic  as  the  pressure  wave  travels,  H  is  the  enthalpy 
and  can  be  determined  from  the  bubble  wall 


In  the  enthalpy  expression,  A  —  0.3001  GPa,  B  =  0.3  GPa, 
n  =  7,  po  =  1000  kg/m3,  and  cq  =  1480  m/s, 

if  we  consider  water  is  incompressible,  the  pressure  field 
can  also  be  given  by  the  following  simplified  expression  [24]: 


P(r)  =  P0  +  —  (p<  -  PO)  + 
r 


Rh(r'  -  Rf,)  ppU2 


r<  2  ('2) 

where  Pc  is  the  bubble  pressure  at  its  radius  Rt,  r  is  the 
distance  from  the  bubble  center,  and  U  is  the  speed  of  the 
bubble  when  the  bubble  expands  to  The  kinetic  energy 
carried  by  the  pressure  wave  can  be  calculated  by  integrating 
the  pressure  pulse  over  the  duration  r'  at  a  given  distance  rs 


Es  = 


PQC0  . 


dl. 


(13) 


C  Power  Dependence 

In  electrical  discharges,  the  energy  delivered  to  the  water 
spark  is  calculable  from  the  measured  voltage  and  current. 
The  generation  of  shock  wave  mostly  occurs  in  the  current 
rise  time  during  which  the  momentum  of  the  expansion  is 
the  largest.  At  later  time,  the  energy  deposition  in  the  plasma 
channel  becomes  much  less,  as  the  majority  of  the  energy  will 
be  absorbed  by  the  external  resistance  in  the  circuit. 

Rather  than  solving  the  differential  equations,  the  peak 
pressure  in  the  plasma  channel  can  also  be  estimated  by  [12] 

p  =  8(J>70*  (14) 

where  P'  is  the  rate  of  power  deposition  in  the  channel  during 
the  resistive  phase  of  the  discharge  and  /  is  the  spark  length. 
The  power  relation  is  simply  given  by  the  product  of  voltage 
and  current 


P(t)  =  I(t)V(t)>  (15) 


HV  Outing  fttwler 


Fig.  2.  Pulsed  power  system  and  experimental  setup.  Two  high  voltage 
capacitors  (I  /*F)  are  connected  in  series  and  switched  by  Tour  MOSFET 
switches  stacking  in  series.  An  exponential-decay  pulse  is  applied  to  the 
electrodes,  which  are  submerged  in  water.  Shock  wave  emits  as  a  result  of 
the  breakdown  of  water  gap. 


HI.  Experimental  Setup 

The  experimental  setup  is  shown  in  Fig.  2,  We  used  two 
needle  electrodes  arranged  in  V-shape  and  the  tips  were  sep¬ 
arated  by  0.2  mm.  The  electrodes  were  mounted  on  a  holder. 
The  two  electrodes  were  touching  the  bottom  of  a  culture 
dish,  which  was  filled  with  tap  water  (resistivity:  —7  kQcm). 
The  pulsed  power  system  used  in  the  experiment  has  two, 
1-pF  capacitors  charged  by  a  dc  voltage  source  (4  kV, 
Ultravolt,  4C  series).  The  capacitors  are  discharged  by  four 
MOSFET  switches  (IXFB38N100Q2)  stacking  in  series.  The 
designing  details  of  such  stacking  circuit  can  be  found  in  [25]. 
In  brief,  the  gate  and  source  terminals  of  the  MOSFET 
switches  are  held  at  the  same  potential  before  closing.  As  soon 
as  the  bottom  switch  is  turned  on  by  a  non  inverted  MOSFET 
driver  (TC4422),  the  second  switch  closes  as  the  gate  is  pos¬ 
itively  biased  through  discharging  the  capacitor  that  connects 
to  the  gate.  In  a  similar  fashion,  the  third  and  fourth  switches 
turn  ON  sequentially.  Because  the  rated  current  of  each 
MOSFET  is  20  A  (average  rating,  the  peak  current  is  higher 
as  shown  later),  several  resistors  connecting  in  series  with  the 
needle  electrodes  were  used  to  limit  the  current.  Four  resistors 
with  a  total  resistance  of  100  £1  were  used.  The  voltage 
across  the  electrodes  was  measured  by  a  high  voltage  probe 
(Tektronix  P6015)  and  the  current  was  measured  by  a  Pearson 
coil  (Pearson,  5046).  A  separate  current  viewing  resistor 
(1.1  Q)  was  also  used  to  measure  the  current  and  was 
connected  in  series  with  the  water  gap.  The  oscilloscope  was 
Tektronix  TDS640A. 

The  laser  shadowgraphy  was  used  to  capture  the  shock 
waves  and  bubble  expansions  after  the  breakdown  of  water 
gap.  A  He-Ne  laser  was  used  to  illuminate  the  electrode 
region  and  a  fast  camera  (4  Pi  cos,  Stanford)  was  used  to 
record  the  images  in  a  time-resolved  fashion.  The  camera  was 
triggered  by  the  oscilloscope  in  response  to  the  current  pulse 
rise  through  the  water  gap.  To  shield  the  electric  interference 
from  the  water  breakdown,  the  pulsed  power  system  and  water 
gap  were  placed  in  a  Faraday  cage,  leaving  the  camera  and 
oscilloscope  outside. 
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Fig.  3.  Voltage  measured  across  the  water  gap  and  the  resulting  current.  The 
peak  voltage  is  approximately  —3,3  kV,  The  water  gap  breaks  down  after 
4S  fi$  at  -2.8  kV  at  0  //s  causing  the  increase  of  current.  The  decay  of  the 
current  is  given  by  the  RC  time  constant,  30  ps. 


Fig.  4,  Voltage  and  current  in  a  shorter  time  scale.  The  current  rise  time  is 
approximately  10  ns, 

IV,  Results 

A.  Voltage  and  Current 

The  voltage  across  the  water  gap  after  the  MOSFETs  are 
turned  ON  is  governed  by  the  capacitor  discharge  through 
the  resistor  (1,28  kQ),  The  peak  voltage  is  approximately 
3.5  kV  (Fig.  3)  but  decreases  according  to  the  time  constant 
(640  /is  and  RC  =  1.28  k£2  x  0,5  /iF),  After  approximately 
50  fi  s,  the  water  gap  breaks  down  so  the  current  starts 
flowing  and  it  decreases  exponentially  at  a  new  time  constant 
(50  /.is  and  RC  —  100  Q  x  0.5  ;;F).  This  is  a  statistical  event 
and  the  breakdown  voltage  is  approximately  2.8  ±  0.5  kV 
with  a  delay  time  of  40  ±  10  ^s. 

Fig,  4  shows  a  voltage  and  current  waveform  in  shorter  time 
scale.  In  general,  the  current  increases  from  0  to  10  ns  as  the 
voltage  drops.  The  voltage  and  current  waveform  allows  us 
to  calculate  the  resistance  of  the  water  gap  and  the  power 
deposition.  In  Fig.  5,  the  water  gap  resistance  drops  from 
>2  k£l  to  approximately  900  Q ,  which  occurs  in  less  than 
2  ns  (Phase  I),  The  second  phase  (Phase  11)  of  the  resistance 
drop  occurs  from  —1.4  to  4  ns,  during  which  the  resistance 


Fig,  5,  Resistance  of  the  water  gap  and  the  corresponding  power  deposited 
in  the  gap.  The  data  were  calculated  from  the  voltage  and  current  in  Fig,  4, 
The  dashed  line  is  the  resistance  fitted  as  the  inverse  square  root  of  lime.  The 
two  dotted  lines  mark  the  phase  transitions. 

becomes  the  same  value  of  the  current-limiting  resistance 
(100  £2).  After  that,  Le„  in  Phase  III,  the  resistance  further 
decreases  to  a  negligible  value  until  >10  ns.  Such  three- 
stage  development  of  water  gap  resistance  agrees  with  the 
resistance  drop  reported  in  [23].  It  is  in  the  second  phase 
that  the  largest  amount  of  power  deposited  in  the  water  gap 
(0-4  ns),  which  according  to  Fig,  5,  has  a  peak  power  30  kW. 
We  note  that  in  the  current  and  voltage  waveforms  (Fig,  4),  tie 
phase  boundaries  are  not  so  clear.  Especially,  the  current  lags 
the  voltage  by  approximately  2  ns,  which  suggests  the  water 
gap  is  somewhat  inductive.  Otherwise,  a  pure  resistive  gap 
should  have  exhibited  synchronous  voltage  and  current.  The 
inaccuracy  caused  by  the  inductive  element  can  therefore  be 
estimated  as  2  ns.  In  other  words,  the  temporal  error  of  Phase  l 
transitioning  to  Phase  II  is  approximately  2  ns.  In  Fig,  5.  the 
boundaries  become  clearer  indicated  by  the  slope  change  of 
the  power  deposition  from  phase  to  phase. 

In  (9),  the  resistance  develops  according  to  the  inverse  of 
square  root  of  time.  We  fitted  this  relation  in  Fig.  5  and  drew 
it  as  the  dashed  line 

Ra(t)  =  100  ./fitfl),  (16) 

The  time  t  is  expressed  in  nanoseconds.  This  relation  works 
well  in  the  second  and  third  phases.  Clearly,  the  resistance 
becomes  infinitely  large  as  time  approaches  zero.  However, 
because  the  expansion  of  water  arc  mainly  occurs  in  the  second 
phase,  we  will  ignore  the  process  prior  to  0  ns,  in  which  only 
small  power  is  deposited  in  the  water  gap  as  the  resistance 
is  still  large  compared  with  the  current-limiting  resistance 
(100  Q).  For  our  calculation  of  (I),  (2),  (8),  and  (9),  the 
starting  time  is  0. 1  ns,  at  which,  the  resistance  of  the  wa:er 
gap  is  3 16  Q  (16).  This  value  is  close  to  the  mean  value  of  the 
measured  values  300  dr  45  £2  at  this  particular  instant  (0,1  ns). 

#,  Temporal  Development  of  Shock  Wave  and  Bubble 

The  temporal  development  of  the  postbreakdown  in  shadow 
graphs  is  shown  in  Fig,  6.  Two  events,  shock  wave  emission 
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Fig,  e.  Laser  shadowgraphs  cf  the  water  gap  after  breakdown.  The  lime 
showr  is  the  delay  from  the  breakdown. 


Time  (ns) 


Fig,  7.  Shock  wave  radius  taken  from  Fig,  6  at  various  delay  times. 


and  bubble  expansion  are  observed.  The  shock  wave  expands 
in  the  spherical  shape,  which  lasts  from  250  ns  to  4.5  p s. 
Witfin  this  time,  there  is  still  current  flow  in  the  water  gap 
between  the  electrodes  (Fig.  3),  The  radius  o:  the  shock  wave 
is  shown  in  Fig.  7.  The  speed  of  the  shock  wave  is  1 .5  km/s, 
measured  by  the  slope  of  the  radius. 

With  the  disappearance  of  shock  wave  in  the  imaging 
window,  the  bubble  grows  and  reaches  its  maximum  radius  at 
M50  /<s.  It  then  shrinks  and  collapses  into  smaller  bubbles. 
The  radius  of  the  bubble  is  shown  in  Fig.  8*  In  the  growing 
phase,  the  radius  i:i  the  expansion  phase  (<150  /is.  Fig.  8) 
can  be  fitted  as  r  =  ata  using  a  graphic  tool  (Sigma  Plot, 
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Fig.  8.  Bubble  radius  taken  from  the  images  in  Fig.  6.  At  each  data  point, 
N  =  5  images  were  included  and  the  radius  is  shown  as  standard  error  bar. 


regression  wizard).  The  resulted  a  is  0.2  and  d  is  0.5 

r  =  0.2ras.  (17) 

The  unit  of  r  is  millimeter  and  /  is  microseconds.  Alterna¬ 
tively,  r  can  be  in  meter  and  /  in  seconds.  The  time  depen¬ 
dence  f05  is  consistent  with  the  experimental  observation  and 
similarity  analysis  [22]. 

The  time  interval  Tc  between  the  bubble  inception  and 
maximum  diameter  is  given  by  [26] 

''max  —  rc/0.915(/Jo/(poo  —  Pv))°  '  (18) 

where  pr  is  the  vapor  pressure  inside  the  bubble  (2330  Pa  at 
20  °C).  In  Fig.  8,  Tc  is  150  ^s,  which  leads  to  rm ax  =  1.6  mm. 
This  radius  is  close  to  2.2  mm,  measured  in  the  experiment 
(Fig,  8), 

C.  Temporal  Development  of  Water  Gap  Pressure 

The  results  shown  in  Fig.  6  summarize  the  postbreakdown 
phenomena  in  the  water  gap,  but  the  pressure  buildup  in 
the  bubble  (or  plasma  channel)  occurs  earlier.  The  largest 
power  deposition  occurs  in  the  breakdown  phase  (Phase  II), 
as  shown  in  Fig.  5.  Equations  (1),  (2),  (8),  and  (9)  allow  us 
to  calculate  the  pressure  in  the  water  gap  during  this  phase. 
Fig.  9  shows  the  pressure  profile  within  80  ns  starting  from  t  = 
OJ  ns.  The  pressure  increases  to  a  peak  pressure  l  GPa  within 
5  ns.  This  is  consistent  with  Fig.  5,  in  which  the  major  power 
is  deposited  in  the  first  5  ns.  The  pressure,  however,  drops 
rather  rapidly  as  the  water  gap  resistance  becomes  negligible 
compared  with  the  current  limiting  resistance,  causing  less 
power  deposition  in  the  gap.  The  accuracy  of  the  pressure 
calculation  is  corroborated  by  the  close  agreement  of  the 
calculated  current  with  the  measured  current  (Fig.  9). 

The  peak  pressure  of  the  water  gap  can  also  be  estimated 
from  (14).  As  shown  in  Fig.  5,  the  peak  power  deposited  in 
the  water  gap  is  30  kW,  which  occurs  in  approximately  2  ns. 
Considering  the  length  of  the  water  gap.  which  is  0.2  mm,  the 
peak  pressure  was  calculated  as  2.2  GPa.  This  value  is  rather 
dose  to  the  I  GPa  obtained  by  solving  the  energy  conservation 
equations  (l),  (2),  (8),  and  (9). 
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Fig,  9.  Pressure  devdopmenl  in  the  water  gap  in  ihe  early  siage  when 
its  resistance  drops  to  the  value  dose  lo  the  current  limiting  resistance.  The 
calculated  and  measured  currents  show  a  dose  agreement. 


Fig.  10.  Pressure  generated  from  (he  bubble  at  its  maximum  pressure  ( I  GP 
and  5  ns  in  Fig.  9)  over  the  distance.  Two  approaches  were  used:  the  red  line 
is  the  result  of  using  ( 10)  and  the  dashed  line  is  the  result  of  using  ( 1 2),  Also 
shown  is  the  electric  Held  calculated  from  the  electrode  by  applying  a  voltage 
of  4  kV  (the  dotted  line). 


D.  Shock  Pressure 

The  shock  pressure  generated  from  the  expansion  of  bubble 
cart  be  calculated  from  (10)  and  (1 1)  or  separately  from  (12). 
We  used  the  pressure  profile  of  the  bubble  wall  from  Fig,  9. 
The  speed  of  the  bubble  wall  was  calculated  from  (17). 

To  compare  the  results  of  these  two  approaches,  we  calcu¬ 
lated  the  pressure  wave  generated  when  the  bubble  pressure 
reaches  the  maximum  (t  —  5  2  ns  and  P  —  I  GPa  in 
Fig,  9).  Fig*  10  shows  the  pressure  field  outside  the  bubble 
extending  to  a  distance  of  7  mm.  The  pressure  drops  from 
I.I  GPa  to  4  MPa  at  7  mm,  whereas  the  incompressible 
approximation  yields  a  pressure  2.5  MPa.  Both  results  are  very 
close,  suggesting  a  simpler  expression  in  (12)  can  replace  the 
more  lengthy  expression  in  (10)  and  (II). 

Since  Fig.  9  shows  the  development  of  the  pressure  in  the 
bubble,  it  allows  us  to  calculate  the  pressure  in  the  time 
domain,  rather  than  single  time  point  shown  in  Fig.  10. 
We  used  (12)  and  read  the  bubble  pressure  from  Fig.  9, 
the  pressure  temporal  profile  outside  the  bubble  at  different 


Fig,  II,  Temporal  development  of  pressure  ai  different  distances  U-6  mm) 
front  (he  electrode. 


distances  can  be  calculated  in  Fig.  I  I .  We  note  that  the 
pressure  is  for  an  observer  at  different  distances  observing  the 
arrival  of  the  pressure  wave  and  the  propagation  delay  is  not 
included.  In  general,  the  pressure  drops  70%  within  400  ns. 
The  trend  of  pressure  decrease  over  distance  approximately 
follows  the  inverse  of  the  distance.  For  example,  at  l  mm 
(870  ns),  the  pressure  is  7  MPa,  but  it  drops  to  I  MPa  at  the 
distance  of  6  mm. 

The  energy  carried  by  the  pressure  wave  was  estimated 
from  (13).  At  the  distance  of  6  mm,  the  energy  integrated 
in  the  time  frame  (<870  ns  in  Fig,  I  I)  is  1  mj,  which  is  the 
same  as  that  calculated  at  the  distance  of  I  mm, 

V.  DISCUSSION 

The  radius  of  the  bubble  expansion  is  well  described  by 
the  power  relation  (17).  The  duration  from  the  expansion 
to  shrinkage  can  also  be  predicted  by  a  well-established 
relation  (18).  Accordingly,  the  resistance  of  the  water  gap  can 
be  described  by  the  inverse  square  root  relation  (16).  Such 
dependence  is  consistent  with  the  observation  in  other  liquids. 

The  pressure  increase  in  the  bubble  coincides  with  the 
increase  of  current,  or  more  directly,  with  the  increase  of 
power  deposition  in  the  water  arc,  which  occurs  in  less  than 
10  ns.  That  points  to  the  importance  of  using  a  low-impedance 
circuit.  However,  due  to  the  restriction  of  rated  current  of 
MOSFET  switches,  such  low-impedance  circuit  will  have  to 
rely  on  higher  current  switches,  such  as  spark  gap  switches. 
In  our  circuit,  the  time  constant  before  breakdown  from  the 
circuit  is  640  /*s,  yet  the  pressure  decreases  70%  within  the 
400  ns.  Moreover,  the  energy  carried  by  the  pressure  wave  is 
only  I  mJ,  which  is  much  less  than  the  total  energy  stored  in 
the  capacitor  4  J.  Therefore,  we  may  decrease  the  capacitance 
value  and  reduce  the  size  of  the  pulse  generator. 

In  estimating  the  shock  pressure,  we  have  used  two 
approaches.  One  is  to  calculate  from  the  method  of  character¬ 
istics  (10),  and  the  other  is  to  calculate  from  the  approximation 
of  incompressible  liquid  (12).  The  drop  of  the  pressure  wave 
over  distance  is  rather  dose  using  these  two  approaches, 
suggesting  a  strong  dependence  on  the  inverse  of  distance. 
On  the  other  hand,  the  peak  pressure  in  the  bubble,  calculated 
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by  (]),  (2),  (8),  and  (9),  is  1  GPa,  in  the  same  range  as  that 
calculated  by  a  pressure-power  relation  (14),  which  is  2,2  GPa, 

Our  motivation  of  this  paper  is  to  produce  shock  waves 
on  the  coverslip.  We  have  used  electric  discharges  and  it 
takes  finite  length  of  time  to  generate  discharge,  approximately 
50  /is,  as  shown  in  Fig,  3,  The  cells  near  the  electrodes  experi¬ 
ence  not  only  shock  waves,  but  also  pulsed  electric  fields.  Side 
effects  such  as  electroporation  could  be  present  in  addition  to 
the  biological  effects  induced  by  shock  waves.  We  therefore 
have  calculated  the  electric  field  dependence  over  distance 
(the  dotted  line  in  Fig.  10)  for  a  maximum  applied  voltage 
of  4  kV,  The  electric  field  drops  rather  rapidly  over  distance. 
At  a  distance  of  2  mm,  the  electric  field  already  drops  to 
1  kV/cm.  However,  the  threshold  for  electroporation,  as  shown 
by  [27 1,  is  greater  than  1  kV/cm  for  50- /is  pulses.  We  therefore 
rule  out  electroporation  for  cells  situated  beyond  2  mm.  For 
other  biological  effects  considered  less  in  severity,  such  as 
action  potential  activation,  a  separate  control  experiment  can 
be  performed.  In  this  case,  the  water  gap  can  be  enlarged 
slightly  so  that  no  breakdown  occurs.  Cells  will  be  only 
exposed  to  electric  fields  so  the  effective  range  of  shock  wave 
impact  can  be  determined. 

Finally,  we  point  out  one  caveat  of  using  ( I ),  (2),  (8).  and  (9) 
to  calculate  the  pressure.  The  fitted  resistance  of  the  water 
gap  (16)  only  works  well  in  the  second  phase  (0-6  ns),  as 
shown  in  Fig.  5.  Beyond  that,  the  actual  resistance  drops  much 
faster  than  predicted  by  (16).  Thus,  the  power  deposited  in 
the  water  gap  should  be  smaller,  which  suggest  the  calculated 
bubble  pressure  should  drop  faster.  Therefore,  the  duration  of 
the  pressure  (70%)  can  be  shorter  than  400  ns. 

VI.  Conclusion 

We  have  used  a  compact  pulsed  power  generator,  which 
consists  of  four  1-kV  MOSFET  switches,  to  produce  shock 
waves  in  biological  media  for  biomedical  research.  Such 
shock  wave  system  allows  cells  to  expose  to  shock  waves 
on  the  microscope,  therefore  providing  a  means  of  real-time 
measurement  of  cells  response.  The  pressure  increase  occurs 
in  the  rising  time  of  the  current  and  therefore  the  circuit 
impedance  needs  to  be  minimized.  In  general,  pressure  waves 
in  the  range  1-10  MPa  with  a  duration  of  sub  microseconds  can 
be  generated  by  this  device.  We  also  would  like  to  point  out  the 
predicted  pressure  needs  to  be  confirmed  by  measurements. 
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Although  cells  often  can  tolerate  oxidative  environments,  abnormal  oxidative  stress  has  been  identified 
in  inflammation,  cardiovascular  and  neurodegenerative  diseases,  and  aging.  The  impact  of  oxidative 
stress  on  the  cellular  biomechanics  is  poorly  understood,  however.  In  this  study,  we  used  C2C12 
myoblasts  to  investigate  the  effect  of  oxidative  stress,  mimicked  by  hydrogen  peroxide  (H2O2),  on  the  cell 
elasticity  (i.e.,  Young's  modulus),  viability,  and  production  of  intracellular  reactive  oxygen  species  (ROS)+ 
To  better  understand  the  mechanisms  underlying  the  impact  of  H202.  we  examined  various  effectors  of 
the  Rho  signaling  pathway,  which  has  been  shown  to  play  a  key  role  in  the  control  of  cell  mechanics. 
H202  decreased  the  cell  stiffness  in  a  dose-dependent  manner,  caused  cell  death,  and  reduced  the  Rh  >A 
expression  that  was  accompanied  by  down-regulation  of  a-actin,  cytos ke I e ton -membrane  linker 
proteins  (ezrin-radixin-moesion  proteins),  and  focal  adhesion.  Modulating  the  Rho  signaling  by  ushg 
a  Rho  activator  partially  restored  the  cell  stiffness,  enhanced  the  cell  viability,  and  decreased  tie 
intracellular  ROS  level  suggesting  a  potential  intervention  strategy  to  maintain  the  cellular  biomecha¬ 
nical  homeostasis  and  rescue  cell  damage  in  the  threat  of  oxidative  stresses. 

&  2014  Elsevier  ltd.  All  rights  reserved. 


1.  Introduction 

Cells  are  subject  to  oxidative  stress  in  many  situations.  Normal 
metabolic  activities  commonly  involve  redox  processes  (Kurbans 
and  Heintz,  2009:  Gutteridge  and  Halliwell.  2010),  and  thus,  cells 
need  to  cope  with  potential  oxidative  stress  regularly.  Pathological 
situations  arise  when  the  oxidative  stress  in  the  cells  and  their 
microenvironments  exceed  the  coping  capacity,  making  the  in  situ 
cells  susceptible  to  oxidative  damages.  Oxidative  stress  has  been 
associated  with  inflammation  (Tiriball  2005),  cardiovascular  and 
neurodegenerative  diseases  (Hebbel  et  ai.f  1990;  Halliwell,  2000; 
Coyle  and  Pattfarcken,  1993),  and  aging  (Sohal  et  aL  2002),  Other 
examples  of  high  oxidative  stress  involve  post-ischemic  reperfusion, 
such  as  in  the  case  of  neural  tissues  immediately  after  ischemic 
stroke  (Aronowski  et  at,  1997),  cardiac  tissues  after  open-heart 
surgery  (Starkopf  et  aL,  1995),  or  subcutaneous  tissues  at  risk  of 
pressure  ulcers  after  prolonged  excessive  epidermal  loading  that 
compromises  tissue  microcirculation  (Mak  et  at,  2010). 

Cells  in  oxidative  environments  often  need  to  continue  to 
sustain  mechanical  loading.  Tendon  fibroblasts  with  tendonitis 
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continue  to  be  stretched  repeatedly  under  inflammatory  condi¬ 
tions.  Card  to  myocytes  need  to  resume  pumping  blood  immedi¬ 
ately  following  an  open-heart  surgery.  In  the  cascade  of  events 
leading  to  deep  tissue  injury  and  pressure  ulcer  under  epidermal 
loading,  the  skin  fibroblasts  and  keratinocytes  as  well  as  tie 
underlying  adipocytes  and  skeletal  myocytes  may  likely  ie 
reloaded  in  a  re  perfusion -induced  oxidative  environment,  if  tie 
repeated  loading  is  applied  soon  after  the  release  of  the  preceding 
Load-induced  ischemia.  Most  studies  on  the  cell  mechanics  ha^e 
typically  been  conducted  with  the  aim  of  examining  the  normal 
cell  mechanics.  Although  pathological  cell  mechanics  have  been 
reported  for  blood  cells  (Finch  et  ahP  1973;  Wang  et  ak,  2002;  Lam 
et  aL.  2007)  and  cancer  cells  (Cross  et  aL,  2007;  li  et  aL.  2Q0S), 
literature  on  how  oxidative  stress  may  affect  cell  mechanics  and  in 
particular  the  load-carrying  capacity  of  cells  is  scant.  The  RhoA, 
among  others,  is  regulated  by  redox  agents  (Aghajanian  et  cL, 
2009;  Nimnual  et  aL,  2003),  and  this  mech  a  notransduction  pati- 
way  is  expected  to  influence  the  cellular  biomechanics.  We  there¬ 
fore  undertook  this  study  to  elucidate  how  an  increase  in  oxidative 
stress  introduced  by  adding  hydrogen  peroxide  to  the  culture 
media  could  affect  the  cellular  mechanics,  cell  viability,  and 
production  of  reactive  oxygen  species  (ROS)  in  mouse  myoblas:s. 
Mechanistically,  the  Rho  signaling  is  known  to  play  a  key  role  in 
the  control  of  cell  mechanics.  However,  the  role  of  various 
effectors  of  the  Rho  signaling  pathway,  including  actin,  vinculin, 
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ezrin-radixin-moesion  (ERM)  linker  proteins,  remains  to  be  elu¬ 
cidated.  In  addition,  we  further  examined  whether  the  upstream 
regulator  (e.g.,  RhoA)  modulated  these  effectors  and  correlated 
with  cell  death  and,  thus,  suggesting  a  potential  application  of 
biomechanically  strengthening  the  cells  under  the  threats  of 
oxidative  challenges. 

2.  Materials  and  method 

2.L  Cell  culture  and  H202  tjeafrnenr 

Mouse-derived  C2C12  myoblasts  (kindly  provided  to  us  by  Prof.  T-C  He)  were 
cultured  in  complete  media  consisting  of  Dulbecco’s  modified  eagle's  medium 
supplemented  with  \0%  fetal  bovine  serum,  1%  antibiotics/antimycotics  (final 
concentration;  penicillin  100  units/ ml,  streptomycin  100  pg/ml  and  amphotericin 
B  0.25  pg/ml).  in  a  humidified  atmosphere  with  5%  CG2  at  37  C.  with  media 
changed  twice  per  week.  Passages  between  3  to  9  were  used  for  all  experiments. 
For  dose-dependent  treatment  of  oxidative  stress,  C2C12  cells  were  treated  for  24  h 
with  prepared  dilutions  of  H202  in  a  cell  culture  medium  resulting  in  the  final 
concentration  range  from  0  to  2  mM.  These  dosages  were  based  on  the  DNA 
fragmentation  associated  with  apoptotic  cell  death  as  previously  reported  (Siu 
et  aL  2009). 

22  Meosurcttienr  on  cell  stiffness  using  AFM 

An  atomic  force  microscope  (AFM;  Nova  scan  Technologies,  Ames,  1A)  was  used  to 
measure  the  biomechanical  properties  of  individual  mouse  myoblasts.  Full  experi¬ 
mental  details  have  been  provided  elsewhere  (Ti  tush  kin  and  Clio,  2007;  Titushkin 
and  Cho,  2011,  Titushkin  et  a!„  2013;  Sun  el  aL  2013).  Briefly,  to  obtain  a  force- 
distance  curve,  an  AFM  cantilever  [0.12  N/m)  with  a  5-pm  radius  microindentor 
descended  toward  the  cell  at^2mm/s  to  minimize  hysteresis  and  until  a  trigger 
force  of  3  nN  was  reached.  The  force-distance  curves  were  then  analyzed  according 
to  the  Hertz  model  (Radmadier,  2002).  The  average  Young  s  modulus  for  each 
experimental  condition  was  calculated  and  subjected  to  crest  at  the  level  of  a =0,05, 

2.3,  Cell  viability 

C2C12  cells  were  placed  in  a  96-well  plate  to  yield  a  density  of  I  x  1G5  cells/ml 
and  incubated  overnight.  The  original  medium  was  then  replaced  by  fresh  medium 
containing  various  concentration  of  H202  (0-2  mM)  with  or  without  Rho  activator 
II  (Cytoskeleton,  Denver,  CO).  Cell  survival  rate  was  determined  using  the  Cell 
Counting  Kit  8  (CCK-8.  Sigma- Aid  rich.  St  Louis.  MO),  which  provides  a  higher 
detection  sensitivity  than,  for  example,  tetrazolium  salts,  such  as  MTT  or  MTS,  A  10-pl 
CCK-8  solution  was  added  to  the  cells  24  h  post -treatment  and  incubated  for  4  h  at 
37  C.  Absorbance  of  each  well  was  quantified  at  450  nm  using  a  plate  reader 
(BioTek  Instruments  Inc..  Winooski,  VT),  The  cell  survival  was  represented  by  the 
optical  absorbance  under  each  experimental  condition, 

2.4.  Superoxide  labeling 

Mitosox  (Life  Technologies.  Grand  Island,  NY)  was  used  to  label  superoxide 
inside  the  cells.  This  novel,  membrane-permeable  fluorcphore  preferentially 
targets  mitochondria.  Upon  oxidation  by  superoxide,  the  dye  fluoresces  with  the 
emission  maximum  at  580  nm.  Myoblasts  were  treated  with  H202  or  Rho  activator 
II  and  then  incubated  with  1  pM  Mitosox  for  30  min  at  37  C  Cells  were  extensively 
washed  and  immediately  mounted  on  the  microscope  stage  for  imaging. 


3,  Results 

Schematically  shown  in  Fig.  I  A.  individual  cell  stiffness  was 
quantitatively  measured  using  AFM-based  microindentation,  A 
typical  force-deformation  curve  (Fig,  IB)  demonstrates  that  as 
the  AFM  probe  slowly  indents  a  surface,  hard  substrate  quickly 
deflects  the  tip  and  causes  a  sharp  rise  in  the  force-deformation 
curve.  In  contrast,  a  gradual  rise  in  the  curve  is  representative  of 
relatively  soft  substrate,  such  as  the  cell.  Microindentation  results 
showed  that  the  stiffness  of  normal  C2C12  cells  is ^3  kPa.  The  cell 
stiffness  was  found  to  decrease  in  response  to  the  H202  treatment 
for  24  h  in  a  concentration-dependent  manner  (Fig,  1C).  At  the 
H202  concentration  as  low  as  0.2  mM,  the  cell  stiffness  was  only 
slightly  reduced  but  not  statistically  significant.  The  cells  became 


mechanically  softer  (p  <  0,05)  after  the  H202  concentrations  were 
increased  to  >  0,5  mM„  presumably  due  to  cytoskeleta!  reorgani¬ 
zation.  We  also  observed  that  at  the  highest  H2O2  concentration 
(2  mM)  used  in  our  experiments,  cells  were  either  broken  into 
debris  or  shrunk  into  much  smaller  spherical  ceils  that  were  too 
easily  detached.  In  such  a  case,  no  AFM  measurements  could  be 
properly  interpreted  (data  not  shown). 

Although  actins  are  primarily  responsible  for  maintaining  the  cell 
stiffness,  focal  adhesion  can  also  contribute  to  the  cellular  biome¬ 
chanics.  We  therefore  fluorescently  visualized  actin  and  vinculin,  a 
key  protein  found  in  the  focal  adhesion,  and  nuclei  in  response  to  the 
H202  treatment.  In  a  control  experiment,  actin  stress  fibers  were 
noticeably  present,  as  expected.  In  addition,  multiple  focal  adhesions 
were  visualized  in  the  control  cells  and  preferentially  located  at  the 
ends  of  actin  stress  fibers  (Fig.  2A).  Consistent  with  the  cell  stiffness 
measurement  shown  in  Fig.  1  Ct  the  treatment  of  mouse  myoblasts 
with  0,2  mM  H20^  did  not  show  discernible  changes  in  the  actin 
organization  or  focal  adhesions  (Fig.  2B),  At  0.5  mM  H202l  the  actin 
fibers  appeared  slightly  weaker  in  fluorescence  intensity  and  shorter 
in  length  (Fig.  2C).  However,  as  the  H202  concentration  increased  to 
1  mM,  the  actin  stress  fibers  were  no  longer  dearly  delineated,  and 
the  focal  adhesions  as  assessed  by  the  vinculin  punctates  were 
diminished  in  numbers  as  well  as  fluorescence  intensity  (Fig,  2D). 
At  the  2-mM  concentration,  cells  could  no  longer  maintain  the 
proper  morphology,  became  smaller  in  size,  and  readily  detached 
from  the  substrate  (Fig.  2E),  A  quantitative  assay  indicates  a  sig¬ 
nificant  reduction  in  the  vinculin  expression  as  a  function  of  the 
increasing  H202  concentration  (Fig.  2F).  These  findings  suggest  that 
oxidative  damage  mimicked  by  H202  disorganizes  the  actin  stress 
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Fig.  1.  Measurement  of  Young  s  modulus.  (A)  Schematic  drawing  showing  an  AFM 
probe  descending  on  the  cell.  (B)  A  typical  force-deformation  curve  is  obtained 
from  each  AFM  measurement.  (C)  Mean  +  SEM  Young's  moduli  of  60  to  80  C2C12 
cells  24  h  post-treatment  of  varying  H202  concentrations  from  3  to  4  independent 
experiments 
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Fig.  2.  Three-color  imaging  of  C2CI2  ce  Is.  Actin  ( red},  vine uSin  (green],  ard  nucleus  ( blue)  were  fluorescently  visualized  n  control  cells  (A)  and  in  response  to  0.2  (B)r  0.5  (CJ, 
1  (D),  and  2  mM  (E)  H20i  treatment  fer  24  h.  Images  were  acquired  using  a  60  x  objective,  A  quantitative  assay  was  applied  to  measure  the  yinculin  expression  and  to 
represent  the  focal  adhesions  (F). 


fibers  and  focal  adhesions  in  myoblasts  in  a  concentration-dependent 
manner,  and  that  these  subcellular  changes  are  accurately  reflected 
by  AFM-based  micro  indentation  measurement  cf  the  ceil  stfffiess. 

As  ot-actin  is  expressed  in  skeletal  muscle,  we  measured  its 
expression  in  C2C12  myoblasts  using  Western  blot.  The  a-actin 
expression  decreased  as  a  function  of  the  H202  concentrator  and 
its  bands  became  almost  invisible  in  response  to  the  1  or  2  mM 
H2G2  (Fig.  3A)  and  appear  to  correlate  with  the  ce  I  stiffness.  To 
further  examine  the  mechanisms  underlying  the  H2G2-  nduced 
oxidative  stress,  we  measured  the  expression  of  RhoA,  and  the 
family  of  phosphorylated  linker  proteins  (p-ERM).  Simila-  to  the 
changes  in  the  a-actin  expression,  both  the  RhoA  and  p-ERM 
expressions  also  decreased  in  a  concentration-dependent  manner. 
At  the  H2O2  concentration  >  TO  mM,  these  mechanosensitive 
molecules  were  barely  detectable.  Such  decreases  ir  the  proteins' 
expressions  may  be  due  to  cell  detachment  and  subsequent 
death.  We  therefore  repeated  Western  blot  experiments  using 
the  cells  remaining  attached  to  the  substrate  only.  Although  not 
as  pronounced,  both  the  RhoA  and  the  a-actin  expressions 
demonstrated  a  similar  trend  of  down-regularion  11  response  to 


the  H2Q2  treatment  (F  g.  3B).  Using  the  Rho  activator  II  in 
increasing  concentrations  (1-2  pg/ml)t  it  is  apparent  that  the 
expressions  of  RhoA  and  a- actin  were  restored  in  the  presence 
of  oxidative  stress  (Fig.  5C). 

Flydrogen  peroxide  was  found  to  decrease  the  cell  viability  in  a 
concentration-dependen:  manner.  The  cell  viability  was  reduced 
by  approximately  2 -fold  when  myoblasts  were  exposed  to  1  mM 
H^O^,  and  only  an  insignificant  fraction  of  cehs  was  found  viable  in 
response  to  2  mM  H2O2  (Fig.  4A).  We  next  tested  the  role  of  RhoA 
involved  m  the  regulation  of  cell  viability.  A  decrease  in  the  cell 
viability  in  response  to  the  H202  treatment  was  found  to  be 
partially  reversible  by  die  RI10  activation,  which  suggests  that 
RhoA  activation  counteracts  the  effect  of  the  H2G2  oxidative  stress 
(Fig.  4A),  At  as  high  as  l  mM  FI2Q2  concentration,  the  Rho 
activation  appeared  to  provide  some  protective  measures  and 
statistically  improved  the  cell  viability.  At  the  2-mM  concentration, 
virtually  all  cells  were  detached  from  the  substrate  and  no 
interpretable  data  were  obtained,  indicating  that  increased  Rho 
activity  was  not  sufficien :  to  rescue  the  cells  from  severe  oxidative 
stress.  The  rescue  of  cells  from  oxidative  damage  by  Rho  activation 
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Fig.  3,  Western  blot  analysis  of  RhoA,  p-ERM  proteins,  and  ra-actin.  (A)  Treatment 
of  cells  with  various  H202  concentrations  down-regulated  the  RhoA.  p-ERM,  and 
ra-actin  in  a  dose-dependent  manner.  (B]  The  protein's  expressions  were  measured 
using  the  cells  remaining  attached  to  the  substrate  only  following  the  HjOj 
treatment,  (C)  Upreguiation  of  RhoA  and  u-actin  using  the  Rho  activator  II  at 
1  ( +  )  and  2  ( +  + )  pg/ml  in  the  presence  of  1  mM  H^02. 


would  be  expected  to  correlate  with  restoration  of  the  cellular 
biomechanical  properties.  Consistent  with  the  role  of  RhoA  as  a 
media  notransducer,  the  Rho  activation  caused  the  cells  to  become 
stiffer  and,  in  response  to  the  oxidative  stress.  Young's  modulus 
resembled  that  found  in  control  cells  (Fig,  4R).  From  a  mechanical 
point  of  view,  the  Rho  activation  restored  the  biomechanical 
properties  of  C2C12  cells  in  the  presence  of  oxidative  stress. 

We  have  recently  reported  the  family  of  tinker  proteins  (ERM 
proteins)  that  physically  couple  the  cytoskeleton  to  plasma  mem¬ 
brane  are  involved  in  the  regulation  of  cellular  biomechanics  (see 
Titushkin  et  al,t  2010b  for  review).  To  test  the  rote  of  the  ERM 
proteins  in  myoblasts,  siRNA  ERM  was  applied  and  their  expressions 
were  transiently  suppressed  (  >  90%  transfection  efficiency).  Incuba- 
tion  of  C2C12  cells  with  H302  for  24  h  could  not  be  carried  out  when 
the  ERM  proteins  were  knocked  down.  Cells  typically  became 
detached  from  the  substrate  and  a  maximum  of  only  5  h  incubation 
was  experimental  feasible,  suggesting  that  the  cytoskeleton-mem- 
brane  coupling  mediated  by  the  ERM  proteins  regulates  cell  adhesion 
and  likely  cellular  biomechanics.  As  expected,  the  ERM  knockdown 
alone  deceased  the  cell  stiffness  by ^35%  (Fig.  5).  To  be  consistent, 
scrambled  siRNA  was  used  as  control  for  5  h  incubation  with  1  mM 
H203.  When  siRNA  ERM  was  applied  prior  to  exposure  to  the  1  mM 
Ft202  treatment,  no  further  decrease  In  the  cell  stiffness  was 
observed,  implying  that  the  H2O2- induced  cytoskeletal  reorganiza¬ 
tion  (e,g„  cell  softening)  may  be  preceded  by  decoupling  the 
cytoskeleton  and  plasma  membrane. 

One  of  the  likely  signaling  pathways  that  mediate  the  effect  of 
H202  may  be  ROS- induced  ROS  release  (RIRR).  The  mitochondrial 
membrane  potential  could  become  unstable,  and  redox  transition  can 
occur  in  response  to  external  insults  that  typically  lead  to  adverse 
consequences  for  mitochondrial  potential  and  cell  viability.  To  validate 
such  mechanism  in  myoblasts,  we  treated  the  cells  first  with  oxidative 
damage  and  then  visualized  the  production  of  superox  ide  using  a  dye 
(Mitosox)  that  fluoresces  upon  oxidation.  In  control  cells,  the  super- 
oxides  typically  formed  small  clusters  around  the  nucleus  {Figs,  6A 


and  D).  In  response  to  the  1.0-  mM  H202  treatment,  the  expression  of 
superoxide  was  increased,  and  the  dye  was  capable  of  penetrating 
into  the  nucleus  (Figs,  6fi  and  E),  providing  evidence  for  RIRR.  When 
cells  were  treated  with  TO  mM  H2Q2  and  Rho  activator  II.  the 
superoxide  expression  and  distribution  were  partially  restored  and 
resembled  those  observed  in  the  control  cells  (Figs,  6C  and  F), 


4.  Discussion 

In  this  study,  we  investigated  the  role  of  Rho  signaling  in 
response  to  oxidative  stress  using  H202  treatment  and  determined 
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Fig.  4.  Cell  viability  and  stiffness.  (A)  Control  cells  we  e  treated  with  varying 
concentrations  but  without  Rho  activation  (open  bars).  The  Rho  activator  II 
together  with  various  HaQ2  concentrations  were  applied  to  the  cells  (black  bars). 
At  the  2-mM  H2Qz  concentration,  the  cell  viability  could  not  be  reversed  by  the  Rho 
activation,  (B)  Changes  in  the  cell  stiffness  appear  to  correlate  with  the  cell  viability. 
*  indicates  p  <  0.05, 
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Fig.  5.  Modulated  biomechanics  of  C2C12  cells  by  ERM  proteins.  Treatment  of  cells 
with  scrambled  siRNA  did  not  significantly  change  the  cells  stiffness,  as  expected. 
However,  siRNA  ERM  treatment  decreased  Voung's  modulus  by  35^  even  without 
applying  HiOj  (open  bars).  Because  treatments  with  both  and  siRNA  ERM  for 
24  h  could  not  keep  the  cells  attached,  the  incubation  time  was  decreased  to  5  h 
only  (solid  bars).  With  scrambled  siRNA  applied,  the  cell  stiffness  in  response  to  a 
shorter  H20j  incubation  time  decreased  by -20%.  No  further  decrease  was 
observed  following  both  siRNA  ERM  treatment  and  a  5-h  incubation  Data 
represent  mean  ±  SEM  of  60  to  80  C2C12  cells  from  3  to  4  independent 
experiments. 
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Fig.  6.  Prtduaiort  of  intracellular  superoxide.  Mitosox  (1  \aM )  was  used  to  visualize  intracellular  superoxide  in  control  cells  {A*  D)  and  in  cells  treated  with  1  mM  alone 
(B,  E)  and  in  the  combination  of  l  mlV  H*0z  and  1  |i gfml  fiho  activator  li  (C  FJ. 


changes  in  the  cellular  biomechanics  and  cell  viability.  In  addition 
to  actin  disassembly  in  C2C12  myoblasts,  we  also  found  that 
oxidative  stress  could  disrupt  focal  adhesion  and  cytoskeleton - 
membrane  association  in  a  dose-dependent  manner,  and  produce 
ROS  (e.g.,  superoxide).  These  mechanical  events  all  contribute  to 
the  altered  cellular  biomechanics.  Moreover,  the  viability  of  C2C12 
myoblasts  was  protected  by  the  Rho  activity  End  correlated 
with  recovery  in  the  cellular  mechanics  following  an  oxidative 
damage 

H202  is  an  oxidant  that  can  induce  cellular  damage  in  a 
concentration-dependent  manner.  It  may  cause  apoptosis  in 
C2C12  cells  (Siu  et  alM  2009).  induce  fragmentation  of  the  mito¬ 
chondria  reticulum  of  C2C12  cells  (Fan  et  al,  2010),  a  nd  produce  ER 
stress  (Pierre  et  al..  2014).  RIRR  is  now  accepted  as  a  common 
median  sm  for  ROS  generation  that  is  regulated  by  the  cell's 
scavenger  system  to  maintain  the  ROS  homeostasis  (for  review, 
see  Zinkevich  and  Gutterman,  2010).  Biomechanically,  it  is  interest¬ 
ing  to  note  that  athero protective  flow  decreases  intracellular  ROS 
level  while  atheroprone  flow  induces  the  opposite  (Dai  eta 1.,  2007). 
As  the  ROS  generation  appears  to  follow  an  a.tered  biomechanics, 
its  Ultra  zellular  level  may  be  modulated  by  changes  in  the  cellular 
biomechanics.  In  our  study,  all  the  effects  attributed  to  H202  are 
regulated  by  or  depend  on  the  Rho  signaling,  and  the  effect  of 
H202  treatment  on  the  three  mechanotransducers  (RhoA,  ERM 
proteins,  and  a-actin)  should  correlate  in  a  concei  ted  manner  (see 
Hg.  3). 

As  RhoA  is  a  member  of  the  GTPase  family  that  has  been  shown 
to  control  actin  stress  fibers  and  focal  adhesions  (Titushkin  and 
Cho,  2C11;  Titushkin  et  al,  2010a.  2010b).  its  activation  may 
stabilize  actins  and  form  stress  fibers  and  subsequently  regulate 
cellular  cytoskeleton  and  morphology.  The  important  role  of 
cytoskeleton  has  been  repeatedly  demonstrated,  particularly  for 
the  load-bearing  cell  types,  such  as  muscle  cells.  It  is  likely 


involved  in  the  regulation  of  intracellular  sgnaling  to  respond 
dynamically  to  a  changing  extracellular  environment.  This  cytos¬ 
keleton  function  has  often  been  attributed,  among  many  others,  to 
the  Rho  family.  Although  complex  interactions  between  the  Rho 
CTPases  and  downstream  target  molecules  have  been  extensively 
studied,  the  current  study  sheds  light  on  the  role  of  RhoA  on 
several  structural  prote  ns  (vincluin.  ERM  proteins,  and  a-actin)* 
These  m  ec  ha  nosens  it  ive  proteins  likely  coordinate  their  responses 
through  the  Rho  signaling  pathway  when  challenged  by  oxidative 
stress,  and  alter  the  cellular  biomechanics. 

The  cell  membrane  not  only  conforms  closely  to  the  cytoske¬ 
leton,  but  it  also  is  physically  associated  with  the  underlying 
cytoskeletal  structures  (Titushkin  and  Cho,  2009).  Evideice  indi¬ 
cates  that  interactions  between  the  cytoskeleion  and  the  lipid 
bilayer,  mediated  by  specific  linkers  (Le„  ERM  proteins),  are  a 
major  factor  determining  the  mechanical  properties  of  the  cell 
membrane,  such  as  membrane  tension.  The  three  closely  related 
ERM  proteins  serve  as  general  crosslinkers  between  the  cell 
membrane  and  actin  filaments  and  organize  the  cytoskeleton  in 
many  cell  types  (Fehon  and  MaClatchey,  201C;  ^eisch  ard  Fehon, 
2011;  Ai  pin  et  al.  2011),  Both  RhoA-deperdent  and  RboA- 
independent  ERM  activation  mechanisms  have  been  suggested. 
Interestingly,  energy  depletion  of  cells  cause!  membrane  separa¬ 
tion  from  cytoskeleton  and  bleb  formation  (T  tushkln  and  Cho, 
2009),  which  is  also  thoight  be  mediated  by  inhibition  of  the  ERM 
linker  proteins.  Changes  in  the  cellular  biomechanics  through 
interfering  with  the  function  of  ERM  proteins  have  been  observed. 
For  example,  our  labora;ory  reported  that  transient  knockdown  of 
ERM  proteins  can  in  hi  jit  or  promote  stem  c^ll  differentiation. 
Osteogenesis  can  be  significantly  delayed  o^  zven  inhibited  by 
down-regulation  of  EFM  proteins  (Titushk.n  and  Cho,  2011), 
whereas  adipogenesis  appears  to  be  facilitatec  (TitushMn  et  al., 
2013). 
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Fig.  7.  Schematic  for  biomechanics]  changes  by  oxidative  stress.  (A)  Endogenous  RhoA  activity  maintains  proper  biomechanical  homeostasis  in  C2C12  cells,  {B}  External 
factors,  such  as  H^Oj,  adversely  affect  the  RhoA  activity,  which,  in  turn,  reorganizes  actin  stress  fibers,  disassembles  Focal  adhesions,  decouples  the  cytoskeleton  from  the  cell 
membrane  by  inhibiting  the  ERM  proteins,  and  increases  the  intracellular  RQS  production.  All  of  the  processes  contribute  to  the  overall  cellular  biomechanics  and  correlate 
with  the  cell  viability.  Activation  of  the  Rho  molecules  is  postulated  to  compete  against  these  effects,  including  a  decrease  in  the  intracellular  ROS  level  and  partially  restores 
the  cellular  biomechanics  and  enhances  cell  viability. 
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Based  on  our  collective  findings,  we  propose  a  working  model 
as  illustrated  in  Fig,  7.  Although  this  model  is  not  intended  to 
represent  a  complete  list  of  possible  multiple  pathways  that  are 
influenced  by  hydrogen  peroxide- induced  oxidative  stress,  the  cell 
viability  and  altered  biomechanical  properties  of  C2C12  cells  may 
be  regulated  through  the  Rho  molecules,  the  focal  adhesions,  and 
the  linker  ERM  proteins.  The  cell  viability  is  reflected  by  and 
correlated  with  ihe  altered  cellular  biomechanics,  in  addition  to 
actin  stress  fibers,  the  overall  cellular  biomechanics  is  determined 
by  focal  adhesions  and  association  between  the  cytoskeleton 
and  the  plasma  membrane  that  is  mediated  by  the  ERM  proteins. 
Hydrogen  peroxide  could  directly  affect  the  focal  adhesion 
and  ERM  proteins,  or  indirectly  affect  them  through  RhoA. 
Although  this  study  is  not  designed  to  delineate  them,  it  is  clear 
that  the  cellular  biomechanics  can  be  altered  by  one  or,  more 
likely,  a  combination  of  effectors  impacting  the  focal  adhesion 
and  ERM  proteins,  which  contribute  to  regulate  the  biomechanical 
homeostasis  in  response  to  oxidative  stress.  Finally,  the  production 
of  intracellular  RGS  through  the  RIRR  mechanism  may  be 
partially  reversed  by  activating  the  Rho  pathways.  It  is  therefore 
envisioned  that  treatments  targeting  the  Rho  signaling  pathway 
could  potentially  be  applied  to  diseases  caused  by  oxidative  stress. 
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Abstract 

Blast  traumatic  brain  injury  (bTBI)  has  now  been  identified  to  associate  with  adverse  health  consequences  among  combat 
veterans.  Post-traumatic  stress  disorder  linked  with  explosive  blasts,  for  example,  may  result  from  such  brain  injury.  The 
fundamental  questions  about  the  nature,  diagnosis,  and  long-term  consequences  of  bTBI  and  causative  relationship  to 
post -traumatic  stress  disorder  remain  elusive,  however.  A  better  understanding  of  brain  tissue  injury  requires  elucidation 
of  potential  mechanisms.  One  such  mechanism  may  be  generation  of  microcavitation  bubbles  in  the  brain  after  an 
explosive  blast  and  their  subsequent  interaction  with  brain  cells.  Using  a  controlled  electrical  discharge  system,  we  have 
successfully  generated  shock  waves  (-^10  MPa)  and  microbubbles  (20-30 /tm)  in  the  cell  culture  of  mouse  astrocytes. 
Detachment  of  astrocytes  from  the  substrate  after  exposure  to  microbubbles  was  observed,  and  it  depended  on  repetitive 
exposures.  Of  the  cells  that  survived  the  initial  assault,  several  subcellular  changes  were  monitored  and  determined  using 
fluorescent  microscopy,  including  cell  viability,  cytoskeletal  reorganization,  changes  in  focal  adhesion,  membrane  per¬ 
meability,  and  potential  onset  of  apoptosis.  While  the  astrocytes  impacted  by  the  shock  wave  only  demonstrated  essen¬ 
tially  unaltered  cellular  behavior,  the  astrocytes  exposed  to  microbubbles  exhibited  significantly  different  responses, 
including  production  of  reactive  oxygen  species  by  collapse  of  microbubbles.  In  the  present  study,  we  characterized  and 
report  for  the  first  time  the  altered  biophysical  and  subcellular  properties  in  astrocytes  in  response  to  exposure  to  the 
combination  of  shock  waves  and  microbubbles. 

Key  words:  astrocytes;  atomic  force  microscope;  fluorescent  microscopy;  microbubbles;  principle  of  symmetry;  shock 
waves 


Introduction 

According  to  the  Centers  for  Disease  Control  and  Prevention, 
approximately  1,7  million  Americans  sustained  traumatic 
brain  injury  (TBI)  between  2002  and  2006.  Mild  TBI,  which  ac¬ 
counts  for  75%  of  annual  TBI  cases,1’2  has  been  linked  to  pro¬ 
gressive  neurodegeneration  and  neurosomatic  damage  leading  to 
cognitive  deficits.  Blast-induced  mild  TBI  now  represents  a  serious 
challenge  to  military  combat  troops/  One  recent  estimate  suggests 
that  approximately  320,000  troops  suffered  mild  TBI  (e.g.,  con¬ 
cussion)  in  Iraq  and  Afghanistan  by  2008/  Definitive  diagnosis  of 
mild  TBI  is  difficult,  because  symptoms  can  often  be  confused  with 
other  medical  conditions/ 3  Primary  TBI  diagnosis  includes  qual¬ 
itative  assessment  of  the  cranium  and  the  use  of  the  Glasgow  Coma 
Scale  (GCS)  scoring  system  that  relies  on  the  clinician's  observa¬ 
tion  of  verbal  and  motor  skills. 

Intracranial  pressure,  CT  scan,  and  MR1  are  also  used  to  assess 
brain  trauma/  Surprisingly,  CT  imaging  and  MRJ  may  not  detect 


mild  TBl-related  damage/  For  example,  CT  scans  are  abnormal  in 
fewer  than  10%  of  mild  TBI  cases,  and  the  diagnostic  capabilities 
of  conventional  MRI  are  also  limited.  It  is  recently  noted  that  dif¬ 
fusion  tensor  imaging  may  be  able  to  better  quantify  and  aid  in 
diagnosis  of  axonal  injury  or  shearing  trauma.  In  addition,  func¬ 
tional  MRI  can  help  to  identify  physiological  and  biochemical 
parameters  associated  with  mild  TBI,  While  imaging  diagnostics 
continues  to  improve,  a  clearer  understanding  of  the  impact  of  blast 
TBI  at  the  cell  and  tissue  level  would  be  crucial  for  the  development 
of  validated  and  unambiguous  diagnosis  and  subsequent  treatment 
strategies. 

Biological  responses  related  to  and  induced  by  explosive  blasts 
can  manifest  themselves  in  brain  tissue  injury.  Energy  from  shock 
waves  caused  by  explosion  (e.g.,  improvised  explosive  devices)  is 
believed  to  transfer  to  the  brain  through  several  modes/  Com¬ 
pression  of  the  skull  plates  is  possible,  but  other  modes,  such  as 
pressurization  of  intravascular  blood  and  stress  from  rapid  severe 
head  movement,  cannot  be  ruled  out.  Blast  may  push  the  blood  into 
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closed  cranial  vaults  and  canals  and  increase  local  pressures* 
Radiofrequency  and  infrared  energies  may  also  result  from  ex¬ 
plosion,  and  the  electromagnetic  radiation  may  interact  with  highly 
electrically  active  organs  such  as  the  brain.'  9  Accelerated  head 
movement  is  interesting  in  that  shearing  may  be  produced  and  lead 
to  axonal  injury. 

More  likely,  energy  transfer  to  the  brain  is  expected  to  be  me¬ 
diated  by  a  combination  of  the  above  modes  and  those  that  arc  yet  to 
be  discovered.  A  recent  report  by  Goeller  and  associates9  showed 
that  in  a  polycarbonate  ellipsoid  brain  model,  blast  overpressure 
ranging  up  to  !70kPa  produced  cavitation  at  the  conlrecoup. 
Moreover,  it  has  been  suggested  that  the  blast -induced  effects  on 
(he  skull  may  generate  microcavitation  bubbles  thnl  subsequently 
collapse  with  high  pressure,  1  Collapse  of  highly  pressurized 
bubbles  can  generate  penetrating  water  jet  (e.g,  microjet)  that  is 
capable  ol  causing  tissue  injury.  Detection  of  these  microbubblcs 
and  potential  impact  on  the  neural  tissues  remain  virtually  un¬ 
known,  however. 

While  no  blast -induced  microbubblcs  have  been  directly  ob¬ 
served  in  the  brain,  cavitation  of  water  has  been  dearly  demon¬ 
strated.  The  pressure  associated  with  collapse  of  the  microbubble 
can  range  from  0.1  to  2G3VtPa,n  suggesting  microcavitation  col¬ 
lapse  may  potentially  be  responsible  for  tissue  damage.  Several 
head  injury  models12-15  also  indicate  that  formation  of  micro- 
cavitation  bubbles  in  response  to  blunt  force  to  the  brain  is  indeed 
feasible  in  experimental  and  computational  studies.  Elucidation  of 
potential  coupling  mechanism  at  the  cell  level  has  been  lacking, 
however.  To  this  end,  we  have  successfully  engineered  an  electrical 
discharge  system  that  can  generate  shock  waves  and  microbubblcs 
of  different  sizes  and  be  mounted  on  a  microscope  stage  for  im¬ 
aging.  This  unique  and  novel  system  allows  us  to  quantitatively 
examine  subcell  ular  responses  using  biophysical  and  molecular 
techniques.  We  characterized  and  now  report  for  the  first  time 
changes  in  the  biophysical  and  biomechanical  properties  of  astro¬ 
cytes  (major  glial  cell  type). 

Methods 

Cell  cultures 

Mouse  astrocytes  were  purchased  from  the  American  Type  Cul¬ 
ture  Collection  (ATCC  Manassas,  VA).  Cell  were  maintained  in  a- 
modified  Eagle  medium  (a-MEM).  10%  fetal  bovine  scrum  (FBS; 
At  la  n  la  Biologieals,  Lawrenceville  GA),  1%  penicillin/streptomydn 
(Invitrogen,  Grand  Island,  NY).  Astrocytes  were  plated  onto  glass 
coverslips  and  cultured  overnight  before  experiment. 

Fluorescence  staining 

Cell  viability  was  determined  using  LIVE/DEAD  viability  kits 
(Life  Technologies)  that  combine  fluorophores  to  yield  two-color 
discrimination  of  the  live  (green)  and  dead  (red)  cell  population.  As 
a  marker  for  focal  adhesion,  vinculin  was  visualized  using  anti¬ 
bodies.  Briefly,  cells  were  seeded  on  the  top  or  bottom  of  the 
chamber,  exposed  to  shock  waves  and/or  microbubblcs,  fixed  with 
3.7%  formaldehyde,  rinsed  with  phosphate  buffered  saline.  Non¬ 
specific  binding  sites  were  blocked  using  a  5%  bovine  serum  al¬ 
bumin  solution  for  1  h  at  room  temperature,  and  then  cells  were 
treated  with  mouse  monoclonal  antibodies  against  vinculin  at  I  :HX) 
dilution.  After  incubation  overnight  at  4 °C,  cells  were  further  in¬ 
cubated  with  FI TC- conjugated  goat  anti-mouse  IgG  at  room  tem¬ 
perature  for  I  h.  To  visualize  F-actin,  cells  were  fixed  in  3.7% 
formaldehyde  and  permeabilized  in  cold  (-20&C)  acetone  for  3  min, 
and  stained  with  rhodamine  phalloldin  ( 1:100  dilution,  Invitrogen) 
for  30  min  at  room  temperature. 


Atomic  force  microscopy 

The  cell’s  elasticity  (i.e.,  Young  modulus  or  sometimes  referred 
It)  as  stiffness)  was  measured  using  an  atomic  force  microscope 
(AFM,  Novascan  Technologies,  Ames,  FA)  mounted  on  an  inverted 
Nikon  microscope.16  To  minimize  the  effect  of  glass  substrate  on 
the  cell  elasticity  measurements,  we  used  an  indentation  depth  up  to 
500  nm  (  ^  10%  of  the  average  cell  height)  for  data  analysis.  About 
100  force- distance  curves  were  generated  using  >30  celts  under 
each  experimental  condition.  The  force- distance  curves  were  col¬ 
lected  and  analyzed  according  to  Hertz  model,17  19  which  relates 
the  loading  force  (F)  with  the  indentation  depth  (/)): 


where  v  is  the  Poisson  ratio,  R  is  the  radius  of  the  spherical  indenter 
<  —5  /im),  E  is  the  local  Young  elastic  modulus.  The  cellular  Poisson 
ratio  was  assumed  to  be  0,5,  which  treats  the  cell  as  an  incom¬ 
pressible  material.  Fitting  the  Hertz  model  to  the  experimental 
force  curve  with  a  standard  least  squares  minimization  algorithm 
yielded  the  local  apparent  elastic  modulus  E, 

Membrane  permeability  and  caspase  signaling 

The  FAM-FLICA*  kit  from  Immunochemistry  (Bloomington, 
MN)  was  used  to  detect  the  membrane  permeability  and  activity  of 
the  caspasc  signaling  in  vitro.  This  kit  detects  multiple  caspascs 
including  caspasc- 1 , 2,  3/7, 6, 9,  and  1 0.  Propidiuni  iodide  f  PI }  was 
used  to  lest  the  membrane  permeability.  After  exposure  to  shock 
waves  and/or  microbubblcs,  both  detached  and  adherent  cells  were 
collected  and  washed  with  the  given  buffer  provided  in  the  kit. 
Cells  were  stained  with  the  FAM -PLICA  caspase  inhibitor  reagent, 
PI,  and  Hocchst  33342.  After  incubation  at  37CC  for  30  min,  the 
samples  were  analyzed  using  a  Synergy  HT  plale  reader  (BioTck, 
Winooski,  VT)  to  measure  the  FAM-FLICA.  PL  and  Hocchst 
33342  signals.  Three  independent  tests  were  repeated  for  each 
condition. 

Superoxide  labeling 

Mitosox  (Life  Technologies)  was  used  to  label  superoxide  inside 
the  cells.  This  novel,  membrane  permeable  fiuorophorc  preferen¬ 
tially  targets  mitochondria.  On  oxidation  by  superoxide,  the  dye 
fluoresces  (emission  maximum  580  nm)  and  is  designed  to  visu¬ 
alize  the  production  of  superoxide  using  fluorescent  microscopy. 
Astrocytes  were  exposed  to  the  combination  o:  shock  waves  and 
microbubbles  and  then  incubated  with  1  mM  Mitosox  for  30  min  ai 
37°C.  Cells  were  extensively  washed  and  mounted  on  the  micro¬ 
scope  stage  for  imaging. 

Exposure  chamber 

The  custom-designed  exposure  apparatus  is  shown  in  Figure  I  A. 
T wo  needle  electrodes  were  separated  by  0.7  mm.  The  pulsed  power 
system  used  in  the  experiment  is  two,  1  ^F  capacitors  charged  by  a 
DC  voltage  source  (4  kV,  800  msec  pulses).  The  capacitors  are  dis¬ 
charged  by  four  MOSFET  switches  (TXFB38N1D0Q2)  slacking  in 
series.  Because  the  rated  current  of  each  MOSFET  is  20 A,  several 
resistors  (a  total  resistance  of  1000)  connecting  in  series  with  the 
needle  electrodes  were  used  to  limit  the  current.  The  voltage  across 
the  electrodes  was  measured  by  a  high  voltage  probe  (Tektronix 
P60I5),  and  the  current  was  measured  by  a  Pearson  coil  (Pearson 
Electronics,  Palo  Alto,  CAh  Details  of  the  electrical  discharge  sys¬ 
tem  have  been  rep  tied  elsewhere.26 

The  volume  of  the  chamber  is  2,5  mL  and  was  intentionally 
designed  to  minimize  temperature  rise  w  ithin  the-  chamber  to  <I°C 
(Fig.  IB).  The  principle  of  symmetry  was  exploited  by  placing  the 
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FIG.  I.  Schematics  of  the  experimental  chamber.  (A)  Photo  of  the  chamber  placed  on  an  upright  Nikon  microscope.  (B)  Side  view  of 
the  chamber  with  specific  dimensions  (not  drawn  to  scale).  Experiments  were  performed  using  mouse  astrocytes  (squares)  seeded  at  the 
lop  or  the  bottom  of  the  chamber. 


two  electrodes  the  middle  of  the  chamber  and  plating  cells  either 
4.7  mm  above  or  below  the  electrodes.  The  rationale  is  that  the 
electric  field  of  short  duration  produced  by  electrical  discharge 
pulses  and  shock  waves  should  propagate  equally  in  the  two  op¬ 
posing  directions  to  the  lop  and  bottom  of  the  chamber  where  cells 
are  seeded.  One  important  difference  is  that  the  microbubbles 
should  rise  to  the  lop  of  the  chamber  only.  While  the  cells  seeded  at 
the  top  and  the  bottom  do  experience  electric  field  (estimated  to  be 
negligible,  M).  28  kV/cm  but  only  Tor  <20(1  //sec)  and  shock  wave 
pressure  (estimated  10MPa)>  the  cells  seeded  at  the  top  only  would 
be  exposed  to  further  interaction  with  the  microbubbles.  The 
membrane  potential  was  monitored  in  real  time  to  ensure  that  the 
short  duration  of  electric  field  exposure  did  not  modify  the  cellular 
responses  (data  available  but  not  shown)* 

Statistical  analysis 

Statistical  analysis  was  performed  at  an  alpha  level  of  0.05.  The 
statistical  significance  of  data  was  evaluated  using  analysis  of  vari¬ 
ance*  All  reported  values  represent  mean  ±  standard  error  of  the  mean. 


Results 

Characterization  of  microbubbles 

Microbubbles  generated  by  the  electrical  discharge  system  were 
visualized  by  differential  interference  contract  (DIG)  microscopy 
and  analyzed  by  a  custom-designed  algorithm.  Some  microbubbles 
were  visualized  as  early  as  40  msec  after  the  electrical  pulses  were 
applied,  suggesting  these  initial  bubbles  traversed  at  an  average 
speed  of  MOO  mm/sec.  The  maximum  number  of  microbubbles 
reaching  the  plane  of  seeded  cells,  however,  was  observed  between 
I  and  2  sec  after  the  stimulation. 

We  developed  a  simple  image  analysis  algorithm  to  identify > 
count,  and  characterize  the  number  and  size  of  the  mict  obubblcs  at 
the  Focal  plane  (Fig.  2A),  A  histogram  was  constructed  to  display 
the  distribution  of  the  number  of  microbubbles  and  their  size  (Fig. 
2B).  In  triplicate  experiments,  the  density  of  microbubbles  was 
observed  to  increase  from  1 0.8/mm 2  (I  pulse)  to  1 8. 6/m  nr  (5 
pulses).  Most  of  the  microbubbles  observed  were  less  than  30 /*m 
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FIG.  2,  Characterization  of  microbubbles.  (A)  Differential  interference  contrast  image  of  microbubbles  captured  at  the  focal  plane 
where  cells  were  seeded.  A  custom-designed  software  was  developed  lo  identify,  count,  and  analyze  the  microbubbles,  (B)  A  histogram 
was  constructed  to  demonstrate  distribution  of  microbubble  number  and  size.  Between  148  and  254  microbubbles  were  captured  in 
triplicate  and  analyzed  for  the  three  different  experimental  conditions*  (C)  Because  the  microbubble  size  is  not  normally  distributed,  the 
Tukey  method  was  applied  to  calculate  statistical  significance.  Data  represent  mean  ±  standard  error  of  the  mean  from  three  independent 
experiments* 
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FIG.  3,  Cell  detachment  and  viability.  Astrocytes  seeded  at  the  top  (A,  C*  F)  or  bottom  (B*  D,  F)  of  the  chamber  were  exposed  to  one 
(A,  B).  three  (C,  D)T  or  five  (E,  F)  shock  waves.  Cell  detachment  was  noticeably  visible  under  some  experimental  conditions  (e,g.,  E; 
five  shock  waves  at  the  top  of  the  chamber).  The  area  of  cell  detachment  (mean  ±  standard  error  of  the  mean)  was  measured  from  three 
to  seven  independent  experiments.  Although  there  was  no  indication  of  cell  detachment  in  the  control  experiments*  (he  areas  of  dark 
regions  were  nonetheless  identified*  and  the  largest  dark  region  was  assumed  to  represent  cell  detachment.  A  lower  magnification 
objective  (20xT  N A/0.5)  was  used  to  capture  large  areas  of  astrocytes.  Color  image  is  available  online  al  www.lieheilpub.com/neu 


radius.  Average  radii  of  the  microbubbles  depended  on  the  number 
of  the  electrical  pulses  applied  (Fig,  2C). 


experiments*  ihc  same  image  analysis  showed  virtually  no  area  of 
cell  detachment. 


Ceil  detachment 


Celt  adhesion  and  biomechanics 


As  expected*  astrocytes  that  were  exposed  to  shcck  waves  and 
microbubbles  became  detached  from  the  substrate,  and  such  ceil 
detachment  was  found  to  depend  on  the  number  of  electrical  dis¬ 
charge  pulses.  Cells  remaining  attached  to  the  substrate  were 
identified  and  visualized  using  a  fluorescent  cell  viabi'ily  assay  that 
fluoresces  live  cells  in  green  and  dead  cells  in  red.  The  cells  ex¬ 
posed  to  one  single  pulse  did  not  show  any  noticeable  detachment 
from  the  substrate  whether  cells  were  plated  at  the  top  or  bottom  of 
the  chamber,  and  the  cell  viability  remained  unchanged  (Fig.  3 A. 
3B),  In  contrast,  a  sequence  of  three  or  five  pulses  exhibited  an 
increasing  extent  of  cel  I  detachment  and  in  the  number  of  dead  cells 
(Fig.  3C-3F). 

It  is  evident  that  the  cell  detachment  and  the  decrease  in  cell 
viability  appeared  to  be  correlated.  This  is  perhaps  best  exemplified 
when  astrocytes  were  exposed  to  five  pulses  (Fig  3E,  3F),  A  large 
area  of  cell  detachment  was  noticeably  visible  at  :hc  top  of  the 
chamber  (Fig.  3E),  and  the  number  of  dead  cells  also  increased*  At 
the  bottom  of  the  chamber  where  tnicrobubblcs  were  prevented  to 
reach,  more  dead  cells  were  found  (Fig,  3F),  One  significant  dif¬ 
ference,  however*  is  that  cells  at  the  bottom  of  the  chamber  es¬ 
sentially  maintained  adherence  to  the  substrate,  anc  no  apparent 
cell  detachment  was  observed.  The  area  of  cell  detachment  was 
analyzed  using  Image!  and  determined  ^32*000  /mr.  hi  control 


Cell  detachment  and  viability  am  strongly  dependent  on  cell 
adhesion  to  the  substrate.  To  examine  altered  cell  adhesion  in  re¬ 
sponse  to  shock  wave  and  microbubble,  focal  adhesions  were  vi¬ 
sualized  by  lluorescently  labeling  vinculin,  and  F-actin  was 
identified  using  phalloidin.  Typical  focal  adhesions  are  known  to 
form  pmictates  that  were  clearly  observed  in  control  cells  (Fig.  4A). 
After  five  pulse  stimulation,  cells  seeded  at  the  bottom  of  the 
chamber  (i.e,  no  microbubbles)  still  maintained  similar  organiza¬ 
tion  of  focal  adhesion  (Fig.  4B).  In  contrast*  cells  seeded  at  the  top 


FIG*  4*  Altered  focal  adhesion.  Fluorescent  images  of  vinculin 
in  astrocytes  seeded  at  the  top  but  without  any  stimulation  (con¬ 
trol;  A)  or  at  the  bottom  (Bj  or  at  the  lop  (C)  of  the  chamber  and 
exposed  five  pulse  simulation.  Image  analysis  showed  that  the 
vinculin  expression  in  cells  in  the  panel  (C)  was  reduced  by 
64,1  ±8.5%.  Data  represent  mean  ±  standard  error  of  the  mean  of 
three  independent  experiments. 
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FIG,  5.  F-actin  organization  and  clastic  modulus  of  astroctyes,  F-actin  was  fluorescenlly  labeled  in  the  control  cells  (A),  in  the  cells 
seeded  at  the  bottom  (B)  or  lop  (Cj  of  the  chamber  (tar  =  2Qmm),  A  schematic  illustrates  the  principle  of  atomic  force  microscope 
(AFM)-based  mieroindcnlalion  (D).  Microindentation  up  to  500 am  depth  is  induced  by  an  AFM  tip  and  the  cantilever  deflection  is 
measured  by  laser  interference  pattern.  Young  modulus  of  astrocyte  seeded  at  the  bottom  of  the  chamber  showed  no  statically  significant 
changes,  while  the  cells  seeded  at  the  top  of  the  chamber  became  mechanically  compromised  (D). 


ol  the  chamber  demonstrated  no  distinctive  punctuates,  but  vinculin 
appeared  to  have  been  down- regulated  and  redistributed  {Fig.  4C). 
Quantitative  analysis  Indies les  lhal  the  vinculin  expression  in  the 
cells  seeded  al  the  top  decreased  by  65%,  Th  s  finding  suggests  that 
mi cro bubbles  caused  astrocytes  to  reorganize  focal  adhesions  and 
therefore  diminished  cell  attachment  lo  the  substrate 

To  visualize  the  microti  lament  organization,  F-ac:in  was  fluor¬ 
escent  ly  labeled  (Fig,  5A-5C),  While  astrocytes  and  other  brain 
cells  arc  not  expected  to  demonstrate  thick  microfilament  stress 
fibers,  the  F-aclin  organization  in  cither  contiol  cells  or  cells  seeded 
at  the  bottom  of  the  chamber  exhibited  F-actin  stress  fibers, 
whereas  in  the  cells  seeded  at  the  lop  of  tfe  chamber,  the  stress 
fibers  appear  to  have  been  diminished.  Loosely  organized  micro¬ 
ti  lament  structure  is  consistent  with  mechanically  softened  cells. 
The  celFs  mechanical  property  was  therefore  measured  using 
AFM,  The  cells  were  microindented  --0,5  /mi  by  an  AFM  tip,16  19 
and  Young  moduli  were  determined  (Fig,  5).  In  control  cells,  the 
elastic  modulus  was  about  0.62  kPa,  approximately  fi  ve- fold  softer 
than,  for  example,  human  mesenchymal  stem  cells  i‘3.2kPai6)  but 
similar  to  other  brain  cell  types  (—0.5  kPa21  *. 

Exposure  of  cells  to  the  shock  wave  only  (e.g.,  at  ihc  bottom  of 
the  chamber)  did  not  significantly  alter  the  biomechanical  proper¬ 
ties,  w  hereas  the  microbubbles  impacting  the  cells  seeded  at  the  top 
of  the  chamber  caused  astrocytes  to  become  mechanically  softer 
(jj<0.05),  suggesting  the  cellular  biomechanical  integrity  may 
have  been  compromised. 

PI  uptake  and  activation  of  caspase  signaling 

Changes  in  the  cellular  biomechanics  may  be  an  indication  that 
the  cell  membrane  permeability  could  have  been  mocified  and 
activation  of  the  a  pop  tot  ic  signaling  pathways  initiated.  To  lest  this 
hypothesis,  we  measured  and  quantified  the  number  of  cells  lhal 
were  permeable  to  PL  Because  PI  is  membrane  :mper meant  and 


excluded  by  viable  cells,  uptake  of  PI  is  indicative  of  structural 
alteration  in  the  cell  membrane.  As  show  n  in  Figure  6A,  only  the 
cells  seeded  at  the  top  but  not  al  the  bottom  of  the  chamber  dem¬ 
onstrated  an  increase  in  the  PI  uptake  in  response  to  microbubbles. 
This  finding  is  consistent  with  an  assay  that  identifies  activation 
of  the  caspase  signaling  pathways  (Fig,  6Bh  a  marker  for  cellular 
damage  leading  to  apoptosis.  It  is  interesting  to  note  that  cellular 
interaction  with  a  lesser  r  umber  of  and  larger  size  microbubbles 
(e.g.,  one  pulse  stimulation)  does  not  appear  to  induce  signifi¬ 
cant  changes  either  in  the  PI  uptake  or  activation  of  the  caspase 
signaling. 

Production  of  superoxioe 

Apoptotic  signaling  can  be  induced  by  oxidative  stress.  To 
measure  the  expression  of  superoxide  in  astrocytes  using  fluores¬ 
cent  microscopy f  cells  were  plated  at  the  top  of  the  chamber  and 
exposed  lo  a  five-pulse  stimulation,  and  immediately  incubated 
with  Mi  to, sox  ( I  mM).  Purctales  of  superoxide  clusters  around  the 
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FIG,  6,  Membrane  permeability  and  apoptosis,  (A)  Uptake  of 
propidum  iodide  (PI)  and  (B)  activation  of  the  caspase  signaling 
pathways  were  quantitatively  determined.  Data  represent 
mean  ±  standard  error  of  the  mean  from  three  independent  ex¬ 
periments.  ’indicates  /><0.05, 
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FIG.  7.  Product  ion  of  intracellular  superoxide.  Milosox  ( I  mM) 
was  used  lo  visualize  intracellular  superoxide  in  control  cells  (A) 
and  in  cells  seeded  at  the  bottom  (B)  or  top  (C)  of  the  chamber 
using  a  20  x,  N A/0, 5  objective.  A  higher  magnification  objective 
(60x,  NA/1.4)  was  used  to  provide  evidence  for  apparent  accu¬ 
mulation  of  the  Milosox  dye  in  the  nuclei  of  astrocytes  that  were 
exposed  to  microbubbles  (D). 


nucleus  were  observed  in  control  ceils  (Fig,  7A),  Similarly,  the 
Milosox  expression  in  astrocytes  plated  at  the  bottom  of  the 
chamber  showed  no  discern  able  differences  when  compared  with 
the  control  cells  (Fig.  7B).  In  contrast,  the  cells  exposed  to  mi- 
crobubbies  not  only  exhibited  a  higher  level  of  su peroxide  but  also 
showed  that  perhaps  Milosox  dye  penetrated  into  the  nucleus  (Fig. 
7C),  For  better  illustration,  a  higher  magnification  objective  was 
used  to  enlarge  astrocytes  in  which  evidence  of  the  dye  accumu¬ 
lation  inside  the  nucleus  was  demonstrated  (Fig.  7D).  Because  the 
dye  is  highly  fluorescent  on  binding  to  nucleic  acid,  it  is  perhaps 
indicative  of  compromised  nuclear  envelope  that  can  alter  the  nu¬ 
clear  trafficking  of  macromolecules. 

Discussion 

A  potential  cause  of  brain  injury  has  been  hypothesized  and 
attributed  to  cerebrospinal  fluid  (CSF)  cavitation.'3""  While  the 
CSF  can  act  as  a  buffer  to  provide  protection  for  the  brain,  the  brain 
and/or  meninges  separate  from  the  skull  on  force  impact,  and  re¬ 
gional  negative  pressures  overcome  the  tensile  strength  of  brain 
tissue.  This  can  create  vacuum  cavity  (Lc„  mienocavitation  bub¬ 
bles),  and  collapse  of  such  cavities  adversely  affects  the  neural 
tissues.33  It  is  still  controversial,  however,  whether  CSF  cavitation 
occurs  during  explosive  blast.  Other  postulates  and  speculations 
have  also  been  proposed,  including  blast -caused  sudden  accelera¬ 
tion  of  the  head  movement  that  decreases  cerebral  intravascular  and 
intracranial  pressure  sufficiently  to  produce  microca vital  ion  bub¬ 
bles  that  collapse  and  may  generate  secondary  shock  waves.10 
While  no  blast- induced  microbubbles  have  been  directly  observed 
in  the  brain,  cavitation  in  water11  or  in  a  surrogate  human  model9 
has  been  clearly  demonstrated. 

The  goal  of  the  current  study  was  therefore  lo  elucidate  the  role 
of  microcavitation  at  the  cell  level  using  in  vitro  astrocyte  culture 
models. 

Because  the  microbubble  generation  system  used  in  this  study 
does  produce  large  but  short-lived  electric  fields  and  shock  waves. 


we  modified  the  exposure  chamber  (Fig.  I B)  to  rule  out  bioeffects 
caused  by  something  other  than  microbubbles.  First,  the  energy 
input  needed  to  produce  shock  waves  may  have  raised  the  tem¬ 
perature  inside  the  chamber.  Simple  calculations  suggest  that  the 
worst  scenario  (K(X)/isee  puise)  in  which  no  heat  diffusion  is  al¬ 
lowed  would  raise  the  chamber  temperature  cPC.  After  a  five- 
pulse  stimulation,  the  temperature  of  the  buffer  h  the  chamber  was 
measured  within  3  sec  and  found  to  be  no  different  than  ambient 
temperature. 

Further,  we  transfected  astrocytes  with  green  fluorescence  pro¬ 
teins  (GFPs),  Because  GFPs  are  sensitive  to  temperature,  they  have 
been  used  lo  calibrate  the  intracellular  temperature.24'2^  We  ob¬ 
served  no  changes  in  the  GFP  signals  in  response  to  any  stimula¬ 
tions  we  applied  (data  available  but  not  shown),  indicating 
temperature  rise  in  the  chamber  and  inside  the  cells  is  negligible. 

Second,  the  results  shown  in  (his  article  were  obtained  using 
astrocytes  that  were  seeded  4.7  mm  vertically  away  from  the 
electrodes.  Numerical  simulation  results  indicate  that  at  this  ver¬ 
tical  distance,  the  electric  field  strength  was  estimated  as  not  in¬ 
significant  but  only  lasts  <200/tsec.  To  experimentally  rule  out  a 
polemia!  effeci  from  the  electric  field,  we  applied  a  fast  reacting 
membrane  potential  dye  (Di -4 -AN  EPFS) 20  and  monitored  its 
fluorescent  intensity.  If  there  were  a  high  electric  field  present  fora 
considerable  period,  the  cell  membrane  should  be  depolarized  and 
the  fluorescent  intensity  would  be  expected  lo  decrease,  Wc  ob¬ 
served  no  such  changes  in  the  membrane  potential  (data  not 
shown),  suggesting  that  the  effects  associated  with  a  short-lived 
electric  field  can  be  ignored. 

Third,  the  shock  wave  pressure  where  the  cells  were  seeded  was 
estimated  at  10  MPa.20  If  the  shock  waves  were  sufficient  to  cause 
cell  detachment,  such  an  effect  should  have  been  observed  both  at 
the  top  and  the  bottom  of  the  chamber.  The  reason  is  that  we 
intentionally  designed  the  exposure  chamber  to  exploit  the  princi¬ 
ple  of  symmetry.  Shock  weaves  produced  at  the  middle  of  the 
chamber  should  propagate  equally  lo  the  top  as  well  as  lo  the 
bottom  of  the  chamber.  However,  cells  exposed  up  to  five  shock 
weaves  only  do  not  show  evidence  of  cell  detachment,  altered  focal 
adhesion,  changes  in  the  biomechanical  properties,  PI  uptake,  or 
activation  of  the  caspasc  signaling  pathways  (Fig,  3F,  4B,  5B,  and 
6,  respectively).  A  logical  conclusion  may  be  drawn  that  transient 
shock  wave  pressure  on  the  order  of  10  MPa  does  not  modulate 
the  astrocytes'  behavior  or  response. 

One  advantage  of  the  chamber  modification  is  to  force  micro- 
bubbles  to  rise  to  the  top  only  and  therefore  delineate  the  bioeffects 
I  hat  can  be  attributed  to  the  microbubbles.  While  it  is  experi  men¬ 
tally  difficult  to  determine  the  lifetime  of  microbubbles,  move¬ 
ments  of  a  few  microbubbles  were  captured  using  high  frame  rates 
of  video  imaging.  Image  analysis  suggests  that  the  microbubbles 
can  move  as  rapidly  as  100 mm/sec.  At  this  velocity,  the  unimpeded 
micro  bubble  movement  would  require  at  least  0,05  sec  lo  rise  lo  the 
lop  of  the  chamber  (  ^  5  mm).  Because  wc  observed  the  maximum 
bubble  density  between  I  to  2  see  after  stimulation,  wc  can  ap¬ 
proximate  (he  lifetime  of  the  microbubbles  between  fractions  of  a 
second  lo  <2  sec. 

Such  an  estimate  is  actually  consistent  with  the  theory  developed 
by  Ferrari  and  coworkers.2  J  The  effects  observed  at  the  top  of  the 
chamber  are  significantly  different.  Unlike  the  cells  seeded  at  the 
bottom  of  the  chamber,  the  cellular  and  subcell  ular  effects  observed 
in  astrocytes  seeded  at  (he  lop  include  noticeable  cell  detachment 
from  the  substrate,  changes  in  focal  adhesion  and  biomechanics, 
and  an  increase  in  the  PI  uptake  and  activation  of  the  caspase 
signaling. 
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Three  possible  explanations  may  be  offered.  First,  we  noted  that 
ihc  density  of  microbubbles  increased  by  70%  with  the  number  of 
pulses  applied  to  generate  them  (Fig.  2B).  and  more  adverse 
bioeffecis  were  observed  with  more  microbubbles.  It  suggests  that 
astrocytes  exposed  to  more  microbubbles  are  prone  to  adverse 
cellular  damage.  Second,  an  alternative  explanation  may  be  due  to 
the  smaller  size  of  microbubble  in  response  to  using  multiple  pulses 
(Fig,  2C),  Tl  is  interesting  to  note  here  that  the  collapse  or  micro¬ 
bubbles  near  cells  is  expected  to  produce  secondary  and  localized 
shock  wave s. 28-30  The  interior  pressure  of  microbubbles  increases 
with  smaller  radii  and  has  been  calculated  using  the  Young-Laplace 
equation,11’31  (AP=  2  s/r)  where  s  is  the  surface  tension  and  r  is  the 
radius.  For  example,  the  pressure  of  I  im  microbubbles  in  water 
can  be  at  least  an  order  of  magnitude  greater  that  the  pressure  in 
J  0  pm  microbubbles. 

Third,  the  collapse  lime  of  microbubbles  also  depends  on  the 
size.  As  demonstrated  by  Takahashi  and  colleagues,30  40  pm  mi¬ 
crobubbles  are  predicted  to  collapse  within  a  few  seconds,  whereas 
10 pm  microbubbles  could  have  a  collapse  time  >100  see,  provid¬ 
ing  a  longer  lifetime  during  which  the  smaller  microbubbles  may 
affect  the  cells.  Collectively,  it  is  feasible  to  postulate  that  the 
combination  of  a  greater  microbubble  density  and  smaller  micro¬ 
bubble  size  with  a  longer  collapse  time  and  higher  interior  pressure 
is  responsible  for  the  adverse  changes  in  astrocytes  documented  in 
this  article. 

Finally,  we  postulate  here  an  alternate  mechanism  that  collapse 
of  microbubbles  produces  hydroxyl  radicals.33’34  The  temperature 
inside  the  collapsing  microbubble  may  be  significantly  elevated  for 
a  short  lime,  causing  pyrolytic  decomposition  and  producing  re¬ 
active  oxygen  species  such  as  hydroxyl  radicals.3  This  potential 
thermochemical  reaction  near  the  cells  may  generate  oxidative 
stress  and  adversely  affect  the  cellular  integrity  and  function.  If 
validated,  the  implication  would  be,  in  addition  to  the  physical 
assault  induced  by  microbubbles,  collapse  of  microbubbles  may 
alter  the  local  chemical  microenvironment  that  further  exacerbates 
adverse  subeel  I  ular  effects,  including  altered  nuclear  transport/16 

Estimating  the  potential  shear  stress  associated  with  the  collapse 
of  microbubbles  remains  a  daunting  challenge.  If  a  microjet  is 
produced  by  collapsing  microbubbles,  the  magnitude  of  induced 
shear  stress  can  be  estimated  by  r  -  p{Vj/h)  where  p  is  dynamic 
viscosity  of  the  medium  (water;  0,89  x  10"4  Pa  s)  and  Vj  is  the  jet 
speed,  and  h  is  the  typical  cross-section  dimension  of  the  jet.37 

In  our  experiments,  most  bubbles  were  measured  in  the  range  of 
20-30  ;nnh  comparable  to  the  cell  size.  The  velocity  of  the  microjets 
in  our  experiment  is  not  available;  however,  for  similar  bubble 
sizes,  it  was  measured  to  be  ^  10m/sec.3K  If  the  jet  diameter  is 
taken  1/10  of  the  bubble  size,  one  can  estimate  the  shear  stress  is 
350  Pa,  As  a  point  of  comparison,  the  critical  shear  stress  for  cell 
detachment  was  reported  to  be  0. 1  kPa34  and  1-3  kPa  lor  membrane 
rupture. Further,  the  recent  estimate  of  collapse  of  bubbles  is 
predicted  to  be  as  high  as  20  MPa. 4  These  estimates  remain  to  be 
experimentally  validated. 

Several  experiments  are  under  way  to  fabricate  microposts  of 
known  mechanical  properties.  It  is  anticipated  that,  on  collapse  of 
microbubbles,  a  shear  stress  will  be  generated  that  can  bend  the 
microposts.  Deflection  of  the  fluorescent  bead-containing  micro¬ 
posts  can  be  monitored,  and  therefore  the  shear  stress  can  be  di¬ 
rectly  measured. 

The  unique  and  novel  exposure  system  we  developed  and 
characterized  for  the  current  study  allows  us  to  quantitatively 
examine  biological  responses.  Using  cell  cultures  of  astrocytes 
(major  glial  cell  type),  we  have  observed  and  quantified  several 


important  cellular  responses  that  are  critical  for  maintaining  the 
cell  viability  and  functionality.  Astrocytes  occupy  between  20% 
and  50%  of  brain  volume  and  exceed  the  number  of  neurons  by 
—  IQ-fold.4041  These  glial  cells  are  also  known  to  regulate  the 
clearance  of  neurotransmi  tiers  and  thereby  prevent  ex- 
cilotoxicity.4'  Results  from  our  experiments  using  such  an  im¬ 
portant  brain  cell  type  might  lead  to  a  better  understanding  of  the 
effect  of  blast- induced  microbubbles. 

Additional  siudies  are  under  way  to  engineer  physiologically 
relevant  three-dimensional  (3D)  collagen  hydrogels  of  co-culiure 
models  using  neurons  and  astrocytes.  A  recently  developed  co- 
culture  model  has  been  applied  to  study  the  common  motor  neuron 
diseases  such  as  amyotrophic  lateral  sclerosis.43  Our  research  ap¬ 
proaches  that  use  state-of-the-art  molecular  and  biomechanical 
techniques  wall  be  applied  to  test  and  validate/refute  the  hypothesis 
that  potential  bioeffecis  in  the  3D  microenvironment  arise  also 
from  localized  membrane  damage,  cytoskeleton  collapse,  necrotic 
cel!  death,  and  activation  of  apoptosis.  Unlike  other  previous 
studies  using  animal  models,  our  shockwave  and  microcavitatton 
generating  system  is  capable  of  repetitive  exposure  to  the  neural 
cells  and  tissues,  A  spatially  and  temporally  integrated  model  may 
be  formulated  in  which  complex  series  of  molecular  interactions  in 
response  to  blast-induced  TB1  may  be  better  explained  by  the 
presence  of  microbubbles. 
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Bubbles  induced  by  blast  waves  or  shocks  are  speculated  to  be  the  major  cause  of  damages  in 
biological  cells  in  mild  traumatic  brain  injuries.  Microbubble  collapse  was  found  to  induce  noticeable 
cell  detachment  from  the  cell  substrate,  changes  in  focal  adhesion  and  biomechanics.  To  better 
understand  the  bubble  mechanism,  we  would  like  to  construct  a  system,  which  allows  us  to  clearly 
differentiate  the  impact  of  bubbles  from  that  of  shocks.  Such  a  generator  needs  to  be  low  profile 
in  order  to  place  under  a  microscope.  A  piezoelectric  transducer  system  was  designed  to  meet 
the  need.  The  system  uses  either  a  flat  or  a  spherical  focusing  piezoelectric  transducer  to  produce 
microbubbles  in  a  cuvette  loaded  with  cell-culture  medium.  The  transducer  is  placed  on  the  side  of 
the  cuvette  with  its  axis  lining  horizontally.  A  cover  slip  is  placed  on  the  top  of  the  cuvette.  The 
impact  of  the  waves  to  the  cells  is  minimized  as  the  cover  slip  is  parallel  to  the  direction  of  the  wave. 
Only  bubbles  from  the  medium  reach  the  cover  slip  and  interact  with  cells.  The  effect  of  bubbles 
therefore  can  be  separated  that  of  pressure  waves.  The  bubbles  collected  on  a  cover  slip  range  in 
size  from  100  jum  to  10  pm  in  radius,  but  the  dominant  size  is  20-30  pm.  ©  2015  A!P  Publishing 
LLC  [ http://dx.doi.org/10. 1 063/1 .4936555] 


I.  INTRODUCTION 

Soldiers  who  experience  the  repeated  impact  of  blast 
waves  may  develop  mild  traumatic  brain  injuries  (TBls),  for 
example,  in  the  form  of  a  concussion.*  While  the  cause  of 
TBI  is  still  not  well  understood,  a  leading  theory  is  that 
microbubbles  are  generated  due  to  the  pressure  impacts, 
which  subsequently  collapse,  causing  microstreams.2  The 
force  created  by  the  microstreams  could  reach  a  level  sufficient 
to  disrupt  cell-cell  connections  and  even  create  membrane 
pores. M  This  theory  has  been  proved  in  astrocytes,5  which 
provides  insights  for  3-D  tissue  culture  or  animal  studies.  In  a 
tissue  subjected  to  pulsed  high- in  tensity  focused  ultrasound 
(pHIFU),  acoustic  cavitation  threshold  varies  broadly  and 
depends  on  the  tissue  composition.6  Experimental  shock 
and  blast  systems  that  can  be  repeatedly  applied  include 
electrically  driven  spark  gaps,  focused  laser,  shock  tubes, 
and  piezotransducers.  Each  approach  has  merits  and  disadvan¬ 
tages.  For  example,  a  laser  focused  in  water  can  create  a  single 
shock  and  bubble  precisely  on  the  micrometer  scale,  therefore 
becoming  highly  useful  for  studying  biological  responses  of 
a  single  cell.  In  the  scenario  of  multiple  bubbles,  however, 
several  lasers  are  needed  and  therefore  increase  the  complexity 
of  the  system.3  Shock  tubes  can  easily  reproduce  the  blast 
waves  that  mimic  the  actual  blast  waves  in  a  bombshell 
detonation,  but  due  to  the  large  scale,  their  use  on  an  in  vitro 
microscope  study  becomes  impossible.  Spark  gaps  are  easy 
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to  construct,  especially  in  water  to  produce  bubbles.  In  the 
system  reported  in  our  previous  work,5  7  the  impact  of  bubbles 
on  astrocytes  was  clearly  demonstrated  and  shock  alone  was 
shown  to  be  ineffective  to  induce  cell  damage.  However,  one 
caveat  for  that  system  is  the  presence  of  an  electric  field  before 
the  breakdown  of  the  spark  gap,  which  can  be  a  source  of 
artifacts  and  may  contaminate  the  biological  results.  Although 
we  have  not  observed  the  effect  directly  caused  by  the  electric 
fields  in  astrocytes,  some  other  sensitive  cell  lines  maybe  more 
susceptible  to  activation  by  the  electric  fields  in  addition  to 
the  bubble  effect,  making  the  interpretation  of  results  difficult. 
It  would  be  necessary  to  have  a  cleaner  system  in  which  the 
bubbles  are  the  only  stimuli  to  the  cells.  In  this  paper,  we  report 
such  a  system  with  the  use  of  piezotransducers  immersed 
in  water  to  produce  microbubbles.  The  piezotransducers  are 
driven  by  RF  signals  with  ultra-low  voltages  (peak  voltage 
^40  V)  and  the  electric  fields  are  no  longer  a  concern.  In  the 
Liquid  environment,  the  negative  pressure  in  the  tensile  phase 
allows  bubble  creation  and  growth.  Because  the  transducer 
is  placed  on  the  side  of  the  cuvette  and  the  cover  slip  on 
the  top  of  the  cuvette,  the  pressure  wave's  impact  with  the 
cells  is  minimized.  All  that  interacts  with  the  cells  are  the 
bubbles  rising  to  the  cover  slip  on  which  biological  cells  are 
seeded.  The  biological  response  to  the  bubbles  can  therefore 
be  observed  provided  that  the  microscope  is  focused  on  the 
cover  slip. 

II.  MICROBUBBLE  GENERATOR 
A.  The  electrical  circuit 

We  adopted  a  commonly  used  ultrasonic  water  atomizer 
circuit  to  drive  the  piezotransducer  at  a  resonance  frequency 
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FIG.  I.  A  schematic  diagram  of  the  electrical  drcuti  of  the  microbubble  gen¬ 
erator,  A  tank  circuit  (Ll ,  L3,  and  C3)  determines  the  oscillation  frequency.  A 
power  BJT  Q1  allows  the  replenishment  of  the  energy  stored  in  the  inductors 
fL!  and  L3)  and  sustains  the  oscillation. 


of  -2,1  MHz,  A  schematic  diagram  of  the  driving  circuit  is 
shown  in  Fig.  1.  The  critical  component  is  a  power  BJT  Ql . 
Two  inductors  in  series  (LI  and  L3)  and  a  capacitor  (C3) 
in  parallel  constitute  the  tank  circuit,  which  determines  the 
frequency  of  the  oscillation.  When  Ql  turns  on  as  a  result  of 
the  switching-on  of  Q2  and  charging  of  C3,  L3  and  L  l  become 
charged.  The  voltage  at  LI  feeds  back  through  C4  to  the  base 
of  Ql ,  calling  for  more  current  tor  a  fast  charging  of  L I ,  The 
oscillation  rises  from  the  discharging  and  charging  between 
C3  and  Ll,  which  in  turn  switches  Ql  on  and  off.  In  short, 
Ll,  L3,  and  C3  form  the  oscillation  and  C5  and  C6  provide 
an  AC  short  path  to  the  transducer.  Meanwhile,  the  energy  is 
replenished  by  the  turning-on  of  Ql . 

The  standard  household  voltage  {120  Vrms)  is  stepped 
down  to  -28  Vrms  to  serve  as  the  power  source  of  the  circuit. 
We  note  that  higher  voltages,  for  example,  100  Vrms,  can 
cause  visible  electrolysis  and  are  therefore  not  suitable  for 
generating  vibrations.  A  relay  (Tyco  CNS- 35-72}  is  connected 
to  the  primary  side  of  the  transformer  to  control  the  "on- 
time**  of  the  circuit,  thereby  controlling  the  on-time  of  the 
mechanical  vibration  of  the  piezotransducer.  The  on-time  in 
this  study  was  limited  to  -2  s. 


B.  Microbubble  chamber  design 

The  size  of  the  microbubbles  is  typically  between  single¬ 
digit  micrometers  to  hundreds  of  micrometers.  When  they  are 
resting  on  a  glass  surface,  a  conventional  upright  microscope 
is  sufficient  to  observe  these  bubbles.  However,  in  order  to 
incorporate  the  chamber  under  a  microscope  for  biological 
studies,  a  few  criteria  have  to  be  met:  (1)  there  has  to  be  a 
light  pathway  so  that  the  cover  slip  can  be  illuminated;  (2) 
the  overall  chamber  size  (in  particular  its  height)  should  be 
small  enough  so  that  it  can  be  placed  on  the  microscope  stage; 
(3)  there  should  be  a  fairly  easy  way  to  place  and  remove  the 
cover  slips;  and  (4)  microbubbles  can  reach  the  cover  slip  with 
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minimal  impact  from  the  mechanical  vibration  generated  by 
the  piezotransducer. 

Acrylic  sheet  was  initially  considered  as  the  material 
for  the  chamber  since  they  are  easy  to  machine.  However, 
considerable  heat  damage  to  the  acrylic  sheet  on  the  air 
side  at  the  line-of-sight  of  the  piezotransducer  was  observed 
after  repeated  operation.  This  is  due  to  a  mismatch  between 
the  acoustic  impedance  of  acrylic  and  the  air  interface:  The 
acoustic  impedance  of  acrylic  is  3.26  x  1 0"6  kg  s“ 1  m“2  8  while 
that  of  air  is  400  kg  s" 1  m“2.  Therefore,  the  transition  of  the  low 
impedance  to  high  impedance  causes  an  increased  magnitude 
of  the  mechanical  wave  at  the  interface  and  excessive  heat, 
leading  to  the  damage  of  the  acrylic  sheet  (melting  point: 

-160  °C). 

We  therefore  chose  glass  (melting  temperature  1400  °C) 
as  the  material  to  construct  the  main  body  of  the  chamber.  The 
piezotransducer  was  arranged  at  90&  to  the  cover  slip  to  avoid 
the  direct  transfer  of  vibration  to  the  cover  slip,  A  separate 
compartment  for  the  piezotransducer  was  constructed  with  a 
circular  opening  (O  ~  1 .50  cm)  to  the  primary  chamber.  Edge 
steps  were  constructed  along  the  top  of  the  chamber  walls 
to  support  a  transparent  top  cover,  where  a  stepped  hole  was 
drilled  for  the  placement  of  the  cover  slip.  A  thin  trench  of 
-2  mm  in  width  connecting  to  the  stepped  hole  was  made 
for  easy  mounting  and  dismounting  of  the  cover  slips  with 
tweezers.  The  pre-seeded  cells  on  the  cover  slip  face  the  water 
bulk  in  the  water  so  the  risen  bubbles  interac:  with  them  as  soon 
as  they  make  contact.  Fig.  2(a)  shows  a  schematic  diagram 
of  the  chamber  design  with  primary  dimensions  annotated. 

Fig,  2(b)  shows  a  picture  of  the  chamber  placed  on  Ihe  stage 
of  an  upright  microscope. 

An  added  benefit  of  an  all-glass  chamber  is  that  one  can 
observe  and  photograph  the  generation  of  microbubbles  from 
the  side  of  the  chamber  while  illuminating  the  chamber  with 
a  strong  light  source.  During  our  experiments,  tap  water  was 
used  to  fill  the  chamber  up  to  the  cover  slip.  No  specific 
degassing  procedure  was  carried  out.  If  should  also  be  noted 
that  bubbles  generated  in  previous  cycles  sometimes  attached 
to  the  surface  of  the  piezotransducer.  Ultrasound  energy  in 
the  following  cycle  is  absorbed  by  these  bubbles,  preventing 
further  microbubble  formation  in  water.  One  can  either  extend 
the  duration  of  the  operation,  i.e„  use  the  ultrasonic  vibration 
of  the  piezotrantlucer  to  “shake  off”  the  attached  bubbles, 
or  to  dean  the  surface  of  the  piezotransducer  before  each 
cycle. 

C.  Piezotransducers 

Two  kinds  of  piezotransducer  (STEMiNC)  were  used  in 
this  study.  The  spherical  transducer  has  a  diameter  of  20  mm 
and  a  focal  point  at  30  mm  from  the  transducer  surface.  It 
is  operated  at  the  resonant  frequency  of  2.1  MHz.  The  idea 
was  to  generate  bubbles  preferably  at  the  focal  plane  of  the 
transducer  so  that  the  maximum  amount  of  bubbles  rises  to 
the  cover  slip.  The  fiat  transducer  has  a  diameter  of  19  mm 
and  a  thickness  of  I  mm,  which  has  a  resonant  frequency  of 
2,07  ±  0.05  MHz,  Both  transducers  were  made  with  similar 
materials  and  operate  in  the  thickness  vibration  mode  with  the 
electro-mechanical  coupling  coefficient,  kP  0.59. 
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FIG.  2,  A  schematic  diagram  of  the  microbubblc  generator  chamber  (a)  and 
a  picture  of  the  chamber  under  the  microscope  (b).  The  electrical  wires  that 
connect  to  the  piezotransducer  are  omitted  here.  A  cover  slip  that  has  cells 
pre-seeded  is  placed  in  the  lop  cover*  leaving  the  side  wilh  cells  facing  the 
chamber  water.  The  bubbles  risen  from  the  glass  chamber  interact  with  the 
cells  as  soon  as  they  reach  the  cover  slip. 

III.  CHARACTERIZATION 

A,  Current  and  voltage  characterization 

The  voltage  supplied  to  the  piezotransducer  was  moni¬ 
tored  via  two  voltage  probes  (Tek  P6139A),  The  current 
behavior  was  monitored  through  a  Pearson  current  monitor 
(Model:  2877).  Signals  were  processed  through  a  digital 
oscilloscope  (Tektronix  TDS3054C). 

Fig,  3  shows  a  typical  voltage  waveform  measured  across 
the  piezotransducer,  as  well  as  the  current  through  the  circuit. 
The  peak-to-peak  voltage  is  measured  to  be  about  70  V  with 
the  magnitude  of  the  negative  polarity  slightly  higher  than  that 
of  the  positive.  The  peak-to-peak  current  is  roughly  5  A,  with 
the  positive  current  slightly  higher  than  that  of  the  negative. 
Although  the  driver  circuit  is  of  LC  resonance  in  nature,  but 
the  phase  difference  of  current  and  voltage  is  not  90°  due 
to  extra  stray  components  of  the  transducer  and  connection 
leads.  The  electrical  power  deposition  onto  the  transducer  was 
evaluated  to  be  about  1 7  W.  Considering  the  electromechanical 
coupling  coefficient  of  0,59,  only  about  10  W  is  channeled  into 
the  mechanical  vibration  of  the  piezotransducer. 


FIG.  3.  Currcni  and  voltage  waveforms  of  ihe  circuit.  The  current  is  not  in 
the  same  phase  as  the  voltage.  The  phase  relationship  is  neither  as  that  in  a 
resistive  circuit  nor  a  typical  resonant  circuit. 

690)  in  an  upright  position,  where  the  piezotransducer  was  ar¬ 
ranged  vertically  in  water  (at  about  6  cm  below  water  surface) 
and  the  fiber  tip  (O  =  150  jim)  was  placed  directly  above  its 
surface  at  various  distances.  The  signal  from  the  hydrophone 
controller  was  processed  through  the  digital  oscilloscope, 

A  pressure  waveform  measured  by  the  fiber  optic 
hydrophone  at  a  distance  of  30  mm  from  the  surface  of  the 
spherical  piezotransducer  is  shown  in  Fig.  4.  The  positive 
pressure  reaches  about  9.5  MPa  while  the  negative  pressure  is 
about  -4  MPa. 

The  pressure  of  the  acoustic  wave  in  water  at  various 
distances  from  the  piezotransducer  was  monitored  for  both 
the  spherical  focusing  transducer  and  the  flat  transducer.  The 
results  are  shown  in  Fig.  4(b).  The  spherical  piezotransducer 
was  found  to  have  a  focal  distance  at  30  mm,  which  agrees 
with  the  datasheet  provided  by  the  manufacturer.  The  pressure 
ranges  from  -2  MPa  to  - 1 0  MPa  between  a  fiber-to- transducer 
distance  of  20  mm  and  40  mm.  The  pressure  produced 
by  the  flat  piezotransducer  is,  however,  considerably  lower: 
-1.0  MPa  at  all  distances  measured  (8-40  mm).  Nevertheless, 
both  piezotransducers  were  able  to  generate  microbubbles  in 
water  effectively. 

C.  Imaging 

Images  of  the  bubbles  inside  the  chamber  were  captured 
via  a  digital  SLR  camera  (Nikon  D30GS)  with  a  macro  lens 
coupled  with  a  2x  teleconverter  (Nikon  TC-20E1II).  The  CCD 
of  the  camera  has  a  pixel  resolution  of  181  pixels/mm.  With 
the  2x  lens,  the  camera  should  in  principle  be  able  to  capture 
microbubbles  of  a  diameter  of  roughly  3  pm  (l  pixel).  Bubbles 
grown  large  enough  to  rise  to  the  surface  of  the  cover  slip  due 
to  buoyancy  were  photographed  by  a  CCD  camera  (Olympus 
DP70)  through  an  upright  microscope  (Olympus  BX51)  with 
a  l  Ox  magnification.  The  size  of  the  bubbles  in  the  photos  was 
analyzed  using  Image J? 

IV,  BUBBLE  GENERATION 


B.  Pressure  analysis 

The  acoustic  pressure  generated  in  water  was  monitored  in 
a  separate  chamber  by  a  fiber  optic  hydrophone  (ON  DA  HFO- 
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A.  Bubbles  in  water 

Microbubbles  nucleate  in  the  presence  of  impurities  or 

pre-existing  bubbles  in  water  and  oscillate  in  phase  with  the 
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ptesMirv  ili&trihulian  of  the  spherical  ami  (he  Hat  transducer  flu 

applied  some  wave,  llic  bubbles  expand  during  rare  tact  ions 
(negative  pressure)  mu)  contract  during  compressions  t  positive 
pressure  h  A  typical  si  ill  image  wuli  miembuhbles  in  water  is 
shown  in  Fig.  5m),  The  image  was  taken  at  an  t-stop  ol  4Jfl 
with  an  exposure  lime  ot  40  ms  and  a  synchronized  flashlight 
A  stainless  steel  needle  of  thickness  190  /mi  was  included  as 
a  reference  to  infer  the  size  of  the  bubbles,  11  to  radius  of  the 
bubbles  captured  in  the  images  ranges  ITom  about  10  //m  to 
about  100  /me  By  amity/ ing  a  series  of  images  taken  at  a  rate 
of  24  Iranies/s.  we  observed  that  when  the  pie/oiransdueei  was 
initially  mined  on.  bubbles  generated  in  water  w ere  extremely 
small  ton  I  \  a  few  pixels  in  the  images  r  However,  more  v  isible. 
la  t  get  bubbles  quickly  appeal',  either  due  uuhe  direct  growth  ol 
smaller  bubbles  or  due  to  the  collision  and  merging  of  multiple 
bubbles.  Bubbles  tikis  also  gel  close  to  each  other  and  form 
bubble  clusters  instead  of  merging  into  one  single  bubble. 

A  few  examples  of  bubble  clusters  are  shown  in  Fig,  3(b), 
Interestingly,  no  bubbles  with  a  radius  greater  than  100  //in 
were  observed.  which  max  indicate  the  tipper  limit  ol  the 
bubble  size  in  the  present  system. 

B.  Bubbles  on  cover  slip 

Bubbles  that  grow  large  enough  rise  to  the  surface  ol  the 
cover  slip  and  were  monitored  under  the  microscope  with  a 
I  Ox  objective  lens.  Before  each  run.  bubbles  on  the  cover 
slip  from  i he  previous  cycle  were  cleaned  with  KutiWiprs® 
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delicate  task  w  ipers.  Figure  6b:  is  a  superimposed  image  ot 
17  images  taken  under  the  microscope  at  the  same  location 
ot  the  cover  slip.  The  bold  ol  v  icw  is  rough  I  >  i  2  x  12  mm 
A  total  of  237  bubbles  were  analyzed  through  imagcJ  The 
relay  control  was  set  to  2  s  on-time  and  the  images  were 
taken  after  a  3  s  delay  to  allow  enough  time  for  the  bubbles 
to  rise  to  the  surface  ol  the  cover  slip,  h  should  be  noted  that 
some  small  bubbles  (single  digit  microns  in  diameter)  do  rise 
to  the  surface  of  the  cover  slip  hut  cad  lapse  rather  quickly 
before  other  bubbles  settle  on  t  ic  cover  slip.  The  histogram 
of  the  bubble  radii  is  shown  in  ig  6(b).  The  radii  distribute 
between  3  and  92  //in,  with  the  highest  percent  frequencies 
appealing  between  H)  and  30  /  m.  Bubbles  with  larger  radii 
appeal  less  frequently.  This  resu  t  is  consistent  with  the  size  ot 
the  bubbles  observed  in  wntei  and  that  report  oil  in  Ref.  T  w  here 
the  rnierobubhles  were  produced  by  electiical  discharges  in 
water, 

Bubbles  of  similar  sizes  mi  y  collapse  at  dillcreni  speeds 
on  the  cover  slip.  We  analyzed  a  series  of  896  images  ol 
two  bubbles  captured  in  a  lime  span  of  60  s.  Thirty  images 
with  the  same  time  spacing  were  chosen  and  are  illustrated  m 
Fig.  .  The  two  bubbles  are  noted  as  bubble  I  and  bubble  2, 
respectively,  for  easy  reference  The  initial  radius  of  bubble  3 
{  M  //no  is  bigger  than  that  ol  bubble  2  1 12.5  //in i.  This  lather 
small  difference  in  radii,  however,  leads  to  a  40k  dilfervuce 
in  i he  bubble  volumes  The  tii  ie  evolution  of  the  radii  and 
i lie  calculated  v  olumes  of  the  hi  hides  are  plotted  in  f  igs  sou 
and  §(b).  The  radius  of  bubble  1  decays  much  faster  than  that 

of  bubble  !  Tire  volumes  of  bubbles  I  and  2  decay  steadily 
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FIG,  6.  A  superimposed  image  of  17  images  taken  at  the  same  location  of 
the  cover  slip  (a)  and  the  histogram  of  a  total  of  237  bubbles  tb).  The  bubbles 
whose  radii  distribute  between  10  and  30  //m  appear  with  tie  highest  percent 
frequencies. 


FIG.  B.  The  time  dependence  of  [he  cross-sectional  area  of  the  bubbles  in 
Fig.  7  (a)  and  the  time  dependence  of  the  radius  of  the  cross-sectional  circle 
(b),  These  results  suggest  bubbles  of  similar  sizes  collapse  at  different  speeds 
on  the  cover  slip. 


with  time  and  can  be  described  with  a  3rd  order  polynomial 
function. 

It  is  interesting  to  note  that  the  smaller  the  bubble  gets, 
the  faster  it  collapses  on  the  cover  slip.  Fig.  9  plots  the  radial 
collapsing  speed  versus  the  radii  of  the  bubbles.  The  radial 
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collapsing  speed  of  bubble  1  is  similar  to  that  of  bubble  2 
within  its  detected  radius  range  (8-14  /mi),  while  bubble  2 
collapses  taster  with  the  decrease  of  the  bubble  size. 

This  phenomenon  has  been  consistently  observed  on 
many  bubbles  captured  on  the  cover  slip.  We  generally 
summarize  our  observations  as  follows;  (l)  All  bubbles 
collapse  on  the  cover  slip  given  enough  time.  The  smaller 
the  bubble  is,  the  faster  it  collapses.  Bubbles  with  a  radius  of 
>  12  /jm  take  more  than  60  s  to  collapse  on  the  cover  slip,  (2) 
Bubbles  with  parting  walls  tend  to  coalesce  and  form  bigger 
bubbles,  which  in  turn  collapse  at  their  own  pace.  (3)  Smaller 
bubbles  may  exist  within  larger  bubbles  given  that  they  arrive 
at  the  cover  slip  first. 

Furthermore,  not  all  bubbles  have  the  same  wall  thickness. 
Some  larger  bubbles  with  thick  walls  seem  to  have  smaller 
bubbles  embedded  in  or  attached  to  their  walls  (shown  as 
voids  on  the  dark  walls  of  bubbles  in  Fig,  10),  These  bright 
spots  are  not  all  in  perfect  circular  shape,  which  mean  that  they 
could  also  be  the  results  of  the  collapse  of  smaller  bubbles  on 
the  wall  of  bigger  bubbles. 

V,  DISCUSSION 


FIG.  7.  Photo  montage  of  i wo  similar  bubbles  collapsing  on  the  cover  slip,  ,  . 

The  initial  radius  of  bubble  l  (14  gm)  is  bigger  than  that  of  bubble  2  ^  continuous  sinusoidal  ultrasonic  wave  was  used  in  our 

(12,5/jm).  system.  In  general  pre-existing  gaseous  impurities  serve  as 
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PIC.  9.  Radial  collapsing  speed  versus  the  radius  of  the  cross-sectional 
circle.  As  the  bubble  becomes  smaller,  the  faster  it  collapses. 


the  nuclei  of  bubble  formation.  During  the  negative  portion 
of  the  pressure  wave,  water  is  pulled  apart  at  the  sites  of 
impurity,  forming  acoustic  microbubbles.  During  the  positive 
portion  of  the  pressure  wave,  these  bubbles  are  compressed, 
albeit  not  completely.  This  process  repeats  in  the  following 
acoustic  cycles,  and  the  bubbles  formed  grow  until  reaching 
a  critical  size  known  as  resonance  size .  Beyond  this  size,  the 
bubbles  either  become  unstable  and  collapse  violently  in  a 
few  cycles  or  oscillate  near  the  resonance  size  for  many  more 
cycles.  Bubbles  may  also  become  fragmented  and  reduce  their 
sizes,  which  then  serve  as  nuclei  of  new  cavitation  cycles.  The 
relationship  of  the  resonance  radius  of  the  bubble  with  the 
frequency  /  is  given  by  Young10  as 


where  y  is  the  specific  heat  ratio  of  the  gas  inside  the  bubble, 
which  in  our  case  is  1 ,33  for  steam  (water). 1 3  is  the  ambient 
liquid  pressure,  which  is  approximately  considered  as  1  atm. 
p  is  the  liquid  density,  which  is  999.97  kg/m3  for  water  and 
tu  is  the  angular  frequency  of  ultrasound,  which  in  our  case  is 


FIG  10,  A  photogram  of  bubbles  captured  on  the  cover  slip  with  different 
bubble  wall  thicknesses  and  voids  inside  the  bubbles. 
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2 nf  =  1 ,32  x  1  ()7  Hz.  The  resonance  radius  is  then  evaluated 
to  be  about  1.52  pm.  The  resulting  diameter  of  the  bubbles 
is  -3  pm,  which  is  beyond  the  detection  range  of  our  CCD 
camera. 

However,  two  pathways  may  lead  to  further  growth  of 
these  small  bubbles  beyond  the  resonance  radius  Rr ,  namely, 
rectified  diffusion  and  bubble  coalescence.  Rectified  diffusion 
involves  an  unequal  mass  transfer  across  the  bubble  interface 
in  the  rarefaction  and  compression  phases  of  the  ultrasound 
wave.  Eller  and  Flynn12  suggested  that  gas  diffusion  into/out 
of  a  bubble  depends  on  the  surface  area  of  the  bubble.  A 
bubble  has  a  larger  surface  area  during  its  expansion  than  its 
compression,  therefore  leading  to  more  gas  diffused  into  than 
out  of  the  bubble.  They  also  suggested  that  the  wall  of  bubble 
thins  during  its  expansion,  making  it  easier  for  gas  to  diffuse 
in.  Bubble  coalescence  involves  multiple  bubbles  coming  into 
contact  with  each  other  and  forming  wall  partitions.  The 
parting  walls  rupture  when  they  become  sufficiently  thin, 
leading  to  the  formation  of  bigger  bubbles,13 

Since  our  configuration  has  a  glass  wall  directly  facing 
the  ultrasonic  wave  produced  by  the  piezotransducer  on  the 
tar  end  of  the  chamber,  reflected  wave  interfere  with  incoming 
waves,  leading  to  nodes  and  anti  nodes  of  standing  waves. 
Microbubbles  formed  at  the  antinodes  that  grow  to  a  radius 
>R,  will  be  pushed  toward  the  nodes  by  primary  Bjerknes 
forces  (acoustic  radiation  forces  on  gas  bubbles) 14  and  become 
“inactive  bubbles'’  {bubbles  that  do  not  implode  violently 
within  the  liquid  near  the  site  of  creation).  These  bubbles 
eventually  become  large  enough  and  rise  to  the  surface  of 
the  cover  slip  due  to  buoyancy  once  the  acoustic  wave  is 
interrupted.  A  more  comprehensive  discussion  of  bubbles  in 
an  acoustic  field  can  be  found  in  the  overview  by  Ashokkumar 
and  co-workers.13 

The  radius  of  the  bubbles  that  arrive  at  the  cover  slip  is 
distributed  between  5  and  92  pm,  and  >50%  of  the  bubbles 
have  a  radius  between  10  and  30  pm.  These  bubbles  gradually 
collapse  on  the  cover  slip  or  merge  to  form  larger  bubbles. 
The  smaller  the  bubble  is,  the  faster  the  collapse  is.  We  believe 
these  bubbles  are  of  similar  properties  as  those  produced  by 
the  electrical  discharges  in  water  as  reported  in  Refs.  5  and  7, 
The  collapsing  time  of  bubbles  with  initial  radius  of -12  ^m 
is  about  60  s,  This  long  collapse  time  nay  be  favored  for 
biological  studies  where  the  observation  of  time  evolution  is 
of  essence. 

We  purposely  designed  the  chamber  in  such  a  way  that 
the  direction  of  the  mechanical  wave  is  perpendicular  to  the 
cover  slip,  so  as  to  separate  the  effect  of  the  pressure  wave 
and  the  collapsing  of  bubbles  on  the  cover  slip  during  future 
biological  experiments.  The  all-transparent  chamber  allows 
easy  access  to  imaging  of  the  microbubbles  in  water  through 
the  side  walls  or  observation  of  the  evert  on  the  cover  slip 
under  the  microscope.  On  a  cover-slip  that  has  only  mono- 
layer  cells,  the  dynamics  of  microbubble  may  differ  from  the 
tissue  which  has  different  surface  properties.  Thus,  the  findings 
of  the  results  may  not  be  directly  translatable  to  tissues,  and 
therefore,  our  setup  may  not  be  applied  to  tissue.  But  one  can 
use  3-D  artificial  matrices  to  produce  physiologically  relevant 
multicellular  structures  on  the  cover  slip  and  allow  studies  to 
be  conducted  closer  towards  actual  tissues, 
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Finally.  2.1  MHz  was  chosen  as  the  driving  frequency 
of  the  piezotransducer  in  our  experiments  for  three  reasons: 
(1)  ultrasound  frequency  at  -20  kHz  is  known  to  generate 
large  bubbles  for  mechanical  shearing  applications  such 
as  emulsification:15  (2)  frequencies  between  100  kHz  and 
1  MHz  are  used  to  generate  smaller  bubbles  for  sonochem  ical 
purposes;1*1  and  (3)  ultrasound  with  a  frequency  above 
1  MHz  has  weaker  cavitation  effects  but  is  commonly 
used  for  medical  and  imaging  purposes.1  By  choosing  a 
frequency  of  2.1  MHz*  we  safely  stay  away  from  “power 
ultrasound,”  yet  still  produce  microbubbles  for  biological 
studies. 


VI,  CONCLUSIONS 

An  all-transparent  chamber  was  designed  to  house  a 
vertically  arranged  piezotransducer*  which  was  driven  by  a  RF 
voltage  at  2. 1  MHz  directly  in  water  to  produce  microbubbles 
for  biological  studies  under  a  microscope.  Bubbles  of  radius 
between  5  and  92  ptm  were  observed  on  the  cover  slip  with 
the  highest  percent  frequency  between  1 0  and  30  pm*  These 
bubbles  collapse  with  time:  The  smaller  the  bubble  is,  the 
faster  it  collapses*  Bubbles  of  around  12  pm  in  radius  take 
about  60  s  to  collapse*  The  fact  that  these  bubbles  take 
longer  to  collapse  may  be  beneficial  to  biological  studies  as 
the  mechanical  waves  are  clearly  separated  in  time  from  the 
impact  of  microbubbles,  allowing  easy  analysis. 
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Pitrpose.  Discomfort  from  light  exposure  leads  to  photophobia,  glare,  and  poor  vision  in  71 
patients  with  congenital  or  trauma-induced  Iris  damage.  Commercial  artificial  iris  lenses  are 
static  in  nature  to  provide  aesthetics  without  restoring  the  natural  iris's  dynamic  response  to 
light.  A  new  photo  responsive  artificial  iris  was  therefore  developed  using  a  photochromic 
material  with  self-adaptive  light  transmission  properties  and  encased  in  a  transparent 
biocompatible  polymer  matrix. 

Mmions.  The  implantable  artificial  iris  w-as  designed  and  engineered  using  Photopia.  a  class 
of  photo-responsive  materials  (termed  naphthopyrans)  embedded  in  polyethylene.  Photopia 
was  reshaped  into  annular  disks  that  were  spin-coated  with  polydime  thy  Esiloxane  (PD  MS)  to 
form  our  artificial  iris  lens  of  controlled  thickness.  ?2 

Results,  Activated  by  UV  and  blue  light  in  approximately  5  seconds  with  complete  reversal  in 
less  than  1  minute,  the  artificial  iris  demonstrates  graded  attenuation  of  up  to  40%  of  visible 
and  60%  of  UV  light.  There  optical  characteristics  are  suitable  to  reversibly  regulate  the 
incident  light  intensity,  hi  vitro  cell  culture  experiments  showed  up  to  60%  cell  death  within 
10  days  of  exposure  to  Photopia,  hut  no  significant  cell  death  observed  when  cultured  with 
the  artificial  iris  with  protective  encapsulation.  Nuclear  magnetic  resonance  spectroscopy 
confirmed  these  results  as  there  was  no  apparent  leakage  of  potentially  Loxic  photochromic 
material  from  the  ophthalmic  device. 

Conclusions.  Our  artificial  iris  lens  mimics  the  functionality  of  the  natural  iris  by  attenuating 
light  intensity  entering  the  eye  with  its  rapid  reversible  change  in  opacity  anti  thus  potent  tally  73 
providing  an  improved  treatment  option  for  patients  with  iris  damage.  ?4 

Keywords:  artificial  iris,  iris  functionality,  photo-response  materials,  UV  activation,  polymer 
NMR  spectroscopy 


The  human  iris  provides  eye  color  and  facial  symmetry  while 
serving  a  central  role  in  maintaining  visual  input  to  the 
retina.  Using  the  sphincter  and  dilator  muscles,  the  natural  iris 
quickly  and  dynamically  alters  the  size  of  the  pupil,  allowing  for 
a  30-fold  change  in  light  transmission  to  the  retina.1  ■-  Excessive 
light  exposure  associated  with  diminished  pupil  constriction 
can  lead  to  damage  of  retinal  photoreceptor  cells  causing  a 
decline  in  vision  over  time,*  Congenital  iris  detects,  such  as 
aniridia  or  iris  coloboma,  as  well  as  trauma-induced  iris  defects, 
are  the  two  main  contributors  to  reduced  functionality  of  the 
iris.1 Albinism,  in  which  patients  lack  the  melanin  pigment 
that  provides  opacity  to  the  natural  iris,  can  also  create  a  larger 
area  through  which  excessive  light  can  enter  the  eye.6  The 
reduced  adaptive  response  of  the  iris  to  light  may  result  in 
photophobia,  glare  and  haloes,  double  vision,  and  decreased 
visual  acuity.7 

Treatment  options  to  alleviate  the  adverse  effects  of  iris 
damage  are  limited  to  wearing  sunglasses,  colored  contact 
lenses,  or  artificial  iris  intraocular  implants.  Although  not 
approved  by  the  Food  and  Drug  Administration  (FDA).  Morcher 
GmbH  (City,  State,  Country),  Ophtec  UV  (City,  State,  Country), 
and  Human  Optics  AG  (City,  State,  Country)  fixed-diameter 
artificial  iris  implants  have  been  used  sparingly  in  the  United 
States  under  the  FDA’S  compassionate  care  clause  K  Morcher  S 


and  Ophtec  s  implants  for  partial  or  complete  aniridia  use  a 
fixed  colored  disc  composed  of  polymethylmethacrylate 
(PMMA)f  a  mechanically  stiff  polymer,  which  solely  blocks 
peripheral  light  entering  the  eye.**'11  In  contrast  to  the  hard 
polymers  used  by  Morcher  and  Ophtec,  1  lumanOptics’  artificial 
iris  lens  is  made  of  a  flexible  silicone  polymer  that  allows  for  a 
smaller  incision  site.12  Including  an  image  of  a  natural  iris 
allows  the  Human  Op  tics  implant  to  be  used  as  a  cosmetic 
option  to  restore  or  change  iris  color.  However,  its  fixed  pupil 
size  and  thus  invariable  light  transmittance  fails  to  dynamically 
reduce  bright  light  entering  the  eye.l31t 

Morcher,  Ophtec,  and  Human  Optics'  artificial  iris  implants 
have  been  shown  to  decrease  glare  and  photophobia  as  well  as 
increase  visual  acuity,  but  posisurgica!  complications  include 
decentralization  due  to  migration  of  the  implant,  inflammation 
caused  by  repealed  contact  with  peripheral  tissues,  and 
anterior  uveitis. 1,-17  Additional  undesirable  effects  include 
secondary  glaucoma,  myopic  shift,  and  a  ghosting  phenomenon 
causing  double  vision. |H,! 9  Furthermore,  despite  their  success 
in  reconstructing  the  shape  and  color  of  the  natural  iris,  all 
three  implants  have  a  set  pupil  diameter  and  amount  of  light 
transmitted,  thus  lacking  the  ability  to  adapt  to  changes  in 
incident  light.  One  important  criterion  for  Improving  the 
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current  treatment  options  would  be  to  provide  the  flexibility  lo 
attenuate  the  incident  tight  dynamically. 

The  development  of  photo-responsive  artificial  iris  implant 
should  mimic  the  functionality  of  the  natural  iris  while 
restoring  iris  color  and  shape.  A  well-designed  artificial  iris 
lens  would  continuously  monitor  and  reversibly  decrease  the 
transmission  of  light  based  (m  the  incident  light  intensity.  For 
example,  an  increased  intensity  of  ambient  light  would  result 
in  increased  activation  of  the  photosensitive  component, 
which  should  absorb  and  bk>ck  transmission  of  a  higher 
percentage  of  the  bright  light,  thereby  protecting  the  sensitive 
photoreceptor  cells  of  the  retina.  To  overcome  a  major 
limitation  of  current  iris  defect  treatment  options,  rapid 
activation  and  reversal  times  of  less  than  1  minute  would  he 
preferred  as  well  as  the  synthetic  and  biocompatible  lens 
materials, 

Photochramic  or  photo-responsive  materials  are  activated 
by  light  to  undergo  a  change  in  molecular  conformation.  This 
change  switches  between  a  colorless  more  transparent  state 
and  colored  opaque  state  associated  with  decreased  light 
transmission,20  Photochromic  materials  have  traditionally  been 
used  for  non  medical  applications,  including  memory  storage 
systems,  semiconductors,  and  light  based  electric  switches.  In 
addition,  these  materials  have  been  proven  safe  and  effective  as 
external  ophthalmic  devices  such  as  Transition  lenses. 
Incorporating  a  photo-responsive  material  encased  within  a 
bioeompattbic  inert  polymer  without  leakage  provides  an 
improved  means  of  fulfilling  the  design  criteria.  In  this  study, 
wc  present  a  novel  artificial  iris  that  combines  photochromic 
and  polymer  components  to  create  an  ophthalmic  device  that 
partially  mimics  the  natural  iris  and  can  have  long-term  implant 
stability  in  the  eye,  which  is  a  major  advance  over  current 
commercial  anificial  irises. 


Materials  and  Methods 

Manufacturing  the  Artificial  Iris 

The  artificial  iris  implant  was  fabricated  using  a  photo- 
responsive  material  encased  within  a  polymer  matrix.  Pellets 
of  Photopia,  which  contain  the  photochromic  material  0%  wi/ 
vol  naphthopyran)  in  polyethylene  (Matsui  International, 
Tokyo,  Japan)  were  melted  and  reshaped  using  a  hot  plate 
and  glass  press  to  create  annular  disks  with  the  following 
dimensions:  1  1-mm  outer  diameter,  i-mm  inner  diameter,  and 
1 50-j.tm  thickness  (Fig.  1A).  Encasement  within  a  polymer 
matrix  of  biocompatible  polydimediylsiloxanc  (PDMS;  Sigma- 
A  Id  rich,  St.  Louis,  MO,  USA)  involved  spin  coating  and  curing  a 
lOOpm  layer  on  each  side  of  the  Photopia  disk  to  form  the 
artificial  iris  (Fig.  1 B),  Spin -coaled  thickness  was  reproducible 
within  5%  error.  Polydime thylsiloxane  is  proven  to  be 
biocompatible  and  has  been  widely  used  in  tissue  engineering. 
An  image  of  the  complete  artificial  iris  construct  in  its  photo- 
activated  and  inactive  states  shows  the  distribution  of  Photopia 
within  PDMS  (Fig.  2).  Selection  and  use  of  Photopia,  which  has 
a  distinctive  color  change  from  clear  to  magenta,  induces 
activation  and  subsequent  light  attenuation. 

Optical  Characterization 

The  optical  functionality  of  the  artificial  iris  in  attenuating 
incident  light  was  determined  with  a  UV/Vis  Spectrophotom¬ 
eter  (Beckman  Coulter,  City,  State,  Country).  The  percentage  of 
light  transmitted  through  the  artificial  iris  with  and  without 
activation  by  UV  light  w-as  measured  over  a  range  of 
wavelengths.  Selling  the  wavelength  scan  to  300  to  750  nm 
provided  attenuation  properties  for  the  two  primary  regions  of 


F i Gear  1,  The  light-responsive  artificial  iris  was  fabricated  in  two 
steps:  (A)  photochromic  Photopia  pellets  were  heated  on  a  hot  plate 
and  reshaped  inm  annular  rings  using  a  glass  press  and  cutter,  and  (R)  a 
100-gm  layer  <»F  PDMS  was  spin-coated  on  cadi  side  of  the  artificial  iris 
to  create  a  disk-shaped  construct  of  13mm  diameter  and  350-pm 
thickness. 


interest:  the  UV  (300-400  nm)  and  visible  light  spectrum  (400- 
700  nm).  To  quantify  the  response  of  the  anificial  iris  to 
changes  in  incident  light  intensity,  neutral  optical  density  filters 
(optical  density  0,  0.3,  and  0.6)  ’were  placed  between  the 
Light  source  and  sample  in  the  spectrophotometer.  The 
wavelength  of  maximal  light  attenuation  was  also  determined 
using  these  scans.  Furthermore,  kinetic  scans  at  that  wave¬ 
length  (approximately  360  nm)  led  to  quantification  of  the 
activation  and  reversal  times  for  the  artificial  iris  by  measuring 
the  time-dependent  percent  transmission  changes.  Image 
quality  after  incorporation  of  the  artificial  iris  into  the  path 
of  light  entering  the  eye  was  then  assessed.  A  glass  slide  etched 


A 

Inactive  -  Low  Intensity  Light 


Fnantr  2.  (A)  Design  schematic  of  artificial  iris:  the  annular  dispersion 
of  the  photochromic  material,  Photopia.  is  shown  in  gray.  The  light 
and  tiark  gray  represent  the  change  in  color  and  increased  opacity  of 
Photopia  on  light  activation.  (B)  Images  of  the  artificial  iris  placed  on 
top  of  the  LUC  logo  show  that  visual  transparency  is  maintained  with 
and  without  activation 
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with  uniformly  spaced  lines  and  a  xit)  microscope  objective 
were  used  to  represent  the  object  and  lens  in  a  classic  image 
formation.  Broadening  for  the  thickness  of  the  lines  (c.g., 
blurring)  after  placement  of  the  artificial  iris  between  the  glass 
slide  and  objective  then  served  as  a  measure  to  quantify  the 
light  scattering  expected  after  implantation. 

In  Vitro  Toxicity  Testing 

Primary  human  corneal  fibroblasts  (HCFs),  passage  14  to  16, 
were  cultured  in  Dulbeceos  modified  Eagle's  medium  (Sigma- 
Aldrich)  supplemented  with  10%  fetal  bovine  serum  (Atlanta 
Biologicals,  City,  State,  Country)  and  1%  antibiotic-antimycotic 
solution  (Sigma-Aid rich).  Human  corneal  fibroblasts  were 
seeded  on  22  X  22-mm  coverslips  at  a  cell  density  of  20,000 
cells  per  coverslip.  Indirect  cytotoxicity  experiments  in  which 
HCFs  were  exposed  to  media  containing  the  engineered 
devices  were  conducted*  In  these  tests,  a  glass  covers  I  ip 
seeded  with  HCFs  was  placed  on  lop  of  a  Photopia  disk 
without  encasement  or  the  artificial  iris  in  a  Petri  dish  idled 
with  media.  A  live-dead  cell  viability  assay  (Invitrogen,  City, 
State.  Country)  was  used  to  stain  cells  at  days  L  5,  and  10. 
Experiments  consisted  of  two  samples  each  of  control, 
Photopia,  and  the  artificial  iris.  Three  fluorescent  images  were 
taken  of  each  sample  using  a  confocal  microscope.  Cell  count 
was  obtained  using  Image!  software  ( ha p://imagej, nih.gov/ij/; 
provided  in  the  public  domain  by  the  National  Institutes  of 
Health,  Bethcsda,  MD,  USA). 

Structural  and  Chemical  Characterization  Using 
1JC  NMR 

Tlit  "c  nuclear  magnetic  resonance  (NMK)  experiments  were 
performed  to  elucidate  the  chemical  composition  and  stability 
of  Photopia  and  the  artificial  iris  in  aqueous  solution. 
Specifically  vve  aimed  to  identify  the  reactive,  and  potentially 
toxic,  side  chains  of  Photopia  in  solution.  Photopia  and 
artificial  iris  samples  immersed  in  PBS  for  3  days  were 
evaluated  using  J3C  NMR  along  with  samples  of  supernatant 
PBS  exposed  lor  1  month  to  Photopia  and  the  artificial  iris  to 
determine  potential  leaching  of  substances  in  solution.  The 
NMR  measurements  were  performed  using  [H  decoupled  15C 
NMR  pulse  program  (Brukcr  -  zgpg;  City,  State,  Country)  at 
room  temperature  on  an  8,5-T  (’ll  frequency  360.13  MHz 
and  l3C  frequency  90,55  MHz)  Broker  Advance  spectrometer 
equipped  with  a  quadruple  nucleus  probe  capable  of 
muhinuclear  NMR  measurements.  A  capillary  filled  with  D2C) 
was  placed  in  the  NMR  tube  to  achieve  magnetic  field 
homogeneity  over  the  sample  volume  by  locking  the  deuteri¬ 
um  signal.  The  solid  Photopia  and  artificial  iris  samples  in  PBS 
were  placed  on  top  of  a  doty  aurum  plug  (14  mm)  to  keep  the 
samples  at  the  center  of  the  radiofrequency  coil.  The  90°  pulse- 
width  was  8.75  gs  for  the  monomer  solutions  and  was  used  as 
is  for  polymer  samples.  The  relaxation  delay  was  set  to  2 
seconds  anti  the  number  of  scans  to  1024,  The  experimental 
parameters  for  supernatant  samples  were  the  following:  90° 
pulse  width  13  5  psK  relaxation  delay  2  seconds,  and 
number  of  scans  3072.  The  dimethy  l  sulfoxide  (DMSO) 
carbon  peak  at  39  5  ppm  w7as  used  for  chemical  shift  reference. 

Results 

Characterization  of  Artificial  Iris 

Activation  of  Photopia  and  attenuation  of  incident  light  were 
first  established  through  observation  of  the  distinct  color 
change  from  clear  to  magenta  in  response  to  UV  exposure  (see 


Fig,  2B),  Subsequent  optical  testing  to  quantify  the  function¬ 
ality  of  the  artificial  iris  demonstrated  rapid  but  not  instanta¬ 
neous  activation  vising  UV  and  blue  light  with  maximal 
activation  achieved  in  less  than  36  seconds  (images  not 
shown).  Similarly,  the  artificial  iris  demonstrated  gradual 
reversal  that  began  immediately  after  removal  of  the  light 
stimulus  with  a  return  to  its  original  state  within  I  minute. 
Wavelength  scans  of  the  artificial  iris  showed  maximal 
attenuation  at  390  nm  with  more  than  60%  of  incident  UV 
light  (300-400  nm)  and  50%  to  65%  of  visible  light  (400-700 
nm)  blocked  (Tig,  3A).  It  was  also  noted  that  the  artificial  iris 
attenuated  a  larger  percentage  of  visible  light  in  the  green  to 
orange  (450-650  nm)  range  when  p reactivated  by  UV  light 
exposure  for  30  seconds  (Fig  3A,  UV-activatcd  artificial  iris). 
11  ic  maximum  attenuation  (i,e,,  maximum  absorbance)  was 
found  to  occur  at  390  nm.  I  bis  coincides  with  the  required 
energy  to  reversibly  break  a  carhon -oxygen  bond  in  Photopia 
and  thus  allowing  optical  modulation  of  transmitted  light. 
Characterization  of  the  artificial  irisTs  response  to  changes  in 
incident  light  intensity  was  achieved  using  neutral  density 
fillers  of  optical  density  0,  0,3,  and  0.6  to  modulate  the 
intensity  of  incident  light  by  0%,  50%h  and  75%,  respectively  It 
was  designed  to  create  varying  ambient  light  intensity 
environments.  The  percentage  of  light  blocked  by  the  artificial 
iris  at  a  given  intensity  was  calculated  byr  taking  the  difference 
between  the  incident  light  on  the  artificial  iris  after  traversing 
the  neutral  density  filler  and  the  percent  of  that  light 
transmitted  through  the  artificial  iris.  As  the  optical  density 
of  the  filters  increases  (i.e.,  the  intensity  of  the  light  on  the 
artificial  iris  decreases),  the  percentage  of  that  incident  light 
blocked  by  the  artificial  iris  also  decreases  from  60%  to  10% 
(Fig,  3R),  hi  other  words,  as  the  incident  light  becomes 
dimmer,  the  artificial  iris  transmits  more  light,  thereby 
attenuating  high-intensity  light  while  ensuring  enough  light 
enters  the  eye  in  the  low  light  intensity  environment. 

Experiments  to  assess  post  implant  image  quality  were 
conducted  using  a  glass  slide  etched  with  uniformly  spaced 
lines  100  pm  apart  and  a  X1Q  microscope  objective  to 
represent  an  object-lenslmage  visual  system.  The  control 
image  of  the  glass  slide  was  acquired  using  differential 
interference  contrast  microscopy  (Tig.  iA).  Placement  of  the 
annular  photochromic  region  of  the  artificial  iris  between  the 
glass  slide  and  the  objective  resulted  in  a  darkened  image 
(overall  pixel  intensity  was  reduced  as  seen  by  a  gray 
background  compared  with  the  control  image)  with  blurred 
lines  (Fig,  4R),  The  degradation  of  image  quality  is  perhaps  best 
characterized  by  estimating  the  full  width  half  maximum 
(F’WHM)  values.  For  example,  the  FWHM  was  calculated  5  pm 
in  the  control  experiment  without  introducing  the  artificial  iris 
in  the  optical  pathway.  This  value  can  be  used  to  represent  the 
resolution  and  appears  to  be  independent  of  the  wavelength 
range  in  this  stud;-.  Insertion  of  the  artificial  iris  increased  the 
FWHM  by  approximately  1 0-fold  to  50  pm.  An  adjustment  of 
the  focal  length  by  approximately  150  pm  seemed  to  reverse 
the  blurring  effect  (FWHM  H  pm),  resulting  in  an  image  of 
sharply  refocused  lines  with  decreased  overall  pixel  intensity 
(Fig.  4Q.  Tite  change  in  resolution  of  the  lines  in  the  three 
images  and  the  corresponding  pixel  intensity  modulation  along 
a  horizontal  line  are  shown  in  Figure  41). 

In  Vitro  Response 

Human  corneal  fibroblasts  with  indirect  contact  to  the  artificial 
iris  or  Photopia  over  a  10-day  period  were  stained  with  green 
caleein  AM  (fluorescence  green)  to  identify  live  cells  and 
ethklium  homodimer  to  highlight  the  dead  cells  (red  fluoresces 
red).  The  resulting  images  indicate  that  the  viability  and 
proliferation  of  cells  exposed  to  the  artificial  iris  were 
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Wavelength  (nm) 

Fiiti'ki:  3-  (A)  light  transmittance  of  the  artificial  iris  in  the  ultraviolet  (300-400  nm)  and  visible  (400-700  nm)  spceiruni.  Frcactlvatlon  by 
exposure  to  365-nm  light  for  30  seconds  (UV  activated)  results  in  reduced  transmission  in  the  green  to  orange  (450-650  nm)  range,  thus 
highlighting  the  additional  light  blocked  by  the  artificial  iris  when  outdoors  or  in  the  presence  of  UV  light.  (II)  The  percentage  of  incident  light 
absorbed  by  the  artificial  ids  increases  from  12%  to  60%  as  the  intensity  of  the  incident  light  on  the  construct  is  increased.  Thus,  the  artificial  iris  has 
higher  light  attenuation  in  bright  light  than  it  does  in  dim  light  settings. 


comparable  to  those  of  the  control  sample  (Fig.  5A).  However* 
cells  exposed  to  Photopia  alone  versus  the  control  sample 
showed  decreased  live  cells  after  S  and  10  days.  Quantification 
of  live  HCFs  front  these  fluorescent  images  confirmed  no 
significant  difference  in  the  number  of  live  cells  between  the 
artificial  iris  and  control  samples  at  day  1  through  10  (Fig,  50), 
In  contrast*  on  days  5  and  10,  the  average  live  cell  count  was 
40%  less  for  HCFs  exposed  to  Photopia  compared  with  the 
control 

Structural  and  Chemical  Characterization  Using 
15C  NMR 

Following  our  in  vitro  ex  peri  merits.  NMR  spectroscopy  was 
used  to  determine  Phot  op  la’s  chemical  structure  and  poten¬ 
tially  reactive  side  groups  and  to  validate  the  stability  of  the 
artificial  iris  in  solution.  The  1 3C  NMR  spectra  of  Photopia  and 
the  artificial  iris  immersed  in  PBS  for  3  days  (Fig*  6A)  identified 
the  mobile  molecular  groups  in  the  samples.  The  assignments 
are  given  in  the  Table.  The  common  peaks  seen  in  these  scans 
belong  to  Photopia.  It  was  noted  that  on  coating  with  POMS, 
an  additional  strong  peak  related  to  CH*  of  PDMS  appears  in 
the  artificial  iris  spectra  and  the  peaks  related  to  Photopia  arc 
shifted  upheld.  A  second  set  of  NMR  experiments  (Fig.  6B)  was 
conducted  on  PBS  exposed  to  either  the  artificial  iris  or 
Photopia  over  1  month  to  ascertain  which  molecules  if  any 
leached  out  of  the  samples  over  time.  For  these  samples*  NMR 
of  the  supernatant  liquid  alone  was  tested.  In  the  case  of  no 
leaching.  NMR  peaks  should  not  be  seen,  as  there  are  no 
carbon-based  substances  in  PUS,  However,  both  the  artificial 


iris  and  Photopia  had  CH-O  residual  monomer  peaks,  perhaps 
coming  from  the  main  polymer  anchor  (c,g.,  CH2  of 
polyethylene).  However,  the  Photopia  sample  alone  showed 
an  additional  peak  at  124  ppm  possibly  coming  from  the 
aromatic  carbon  group  in  the  naphthopyran  (Fig,  6(1),  which 
might  be  responsible  for  the  cell  death  in  Photopia  (see  Fig. 
5A)>  This  potential  leaching  is  effectively  blocked  in  the 
artificial  iris,  likely  due  to  PDMS  coating  and  therefore 
improving  the  overall  cell  survival 

Discussion 

The  light-responsive  artificial  iris  was  designed  as  a  biocom- 
paiihlc  ophthalmic  implant  with  self-adjusting  light  transmis¬ 
sion  properties  to  partially  mimic  the  functionality  of  the 
natural  iris.  Dynamic  light  attenuation  was  achieved  through 
incorporation  of  a  photo-responsive  material  (Photopia.  5% 
naphthopyran)  that  undergoes  a  reversible  molecular  confor¬ 
mational  change  induced  by  UVor  blue  light.-1  Photochemical 
cleavage  of  the  carbon-oxygen  bond  within  the  pyran  ring  by 
incident  light  converts  a  naphthopyran  to  its  colored  opaque 
ring-open  state,  while  reversal  to  its  original  transparent  ring- 
closed  form  occurs  spontaneously  with  removal  of  the  UV  light 
stimulus  *2  This  class  of  photochromic  materials  is  particularly 
suited  for  ophthalmologic  applications  due  to  its  efficient 
photo-response,  large  light  absorbance  when  activated,  rapid 
reversal,  and  Fatigue  resistance.23 

Encasing  an  annular  disk  of  Photopia  within  PDMS  (Figs.  I, 
2)  to  form  the  artificial  iris  implant  met  our  design  objectives 
by  restoring  iris  shape  and  providing  significant  light  auenu- 
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Fk;i.sk  4.  (A)  Control.  Image  of  glass  slide  with  black  lines.  (B)  AltiAdal  iris.  Image  through  peripheral  phoCoefafomic  region  of  artificial  iri^  (Al): 

it  gin  same  ring  {increased  thickness  of  blurred  black  lines')  and  decreased  overall  light  transmission  (gray  background)  (O  Artificial  iris  refocused, 
image  through  peripheral  artificial  iris  after  focal  length  increased  by  approximately  150  jam:  visual  acuity  restored.  (D)  Pixel  intensity  versus 
position  along  purple  lines  in  (A)  to  (C),  Image  resolution  denoted  by  the  full  Width  half  maximum  of  control  (A)  and  artificial  iris  refocused  (C)  is 
greater  than  image  formed  through  artificial  iris  alone  (B).  Reduced  light  transmission  (decreased  baseline  intensity)  for  images  through  AI  (B  and  C) 
compared  with  control  (A)  was  also  noted. 


ation  while  maintaining  visual  acuity  through  the  clear  center 
where  light  passes  unimpeded  for  sight.  Through  selection  of  a 
photo-responsive  material  with  activation  and  reversal  limes  in 
less  than  I  minute,  we  have  overcome  the  primary  limitation  of 
other  iris  defect  treatment  options  such  as  Transition  lenses, 
which  have  prolonged  reversal  limes  in  excess  of  10  minutes 
when  indiKirs.  Similarly,  tile  use  of  this  particular  self-adaptive 
Light-sensitive  material  resulted  in  a  Jess  bulky  alternative  to  the 
circuits  found  in  electric  artificial  irises21  and  improved  on 
current  commercial  artificial  irises,  such  as  Mum  an  Optics  AG 
through  its  light-attenuating  profile.21  Encasing  Photo pia  in 
hydrophobic  PDMS  was  central  to  preventing  cell  migration 


and  surface  attachment,  which  could  affect  light  transmission 
through  our  artificial  iris. 

The  human  eye  relies  on  the  pupillary  light  reflex,  the 
constriction  and  recovery  of  the  iris  in  response  to  light 
stimulus,  to  modulate  the  light  flux  into  the  eye/*'  Along  with 
controlling  retina!  illumination,  a  light-responsive  pupil  can 
change  the  depth  of  focus  and,  with  a  smaller  pupil  size, 
reduce  optical  aberrations,  light  diffraction,  and  visual  effects 
such  as  glare/7  2”  Our  goal  in  engineering  a  new  artificial  iris 
was  to  restore  these  natural  iris  properties  to  patients  with 
congenital  or  trauma-induced  iris  defects.  By  blocking  up  to 
60%  of  incident  light  across  the  UV  and  visible  light  spectra, 


P*y  1  D*y  5  ISty  IP 


Figure  5.  Human  comeal  fibroblasts  were  culm  red  with  cither  the  annular  Plmtnpia  disks  or  the  complete  artificial  iris  for  up  to  10  days,  (A)  Cell 
viability  staining  {live  cells  —  green,  dead  cells  —  tvif)  show  similar  morphology  and  proliferation  rate  between  control  and  artificial  iris-exposed 
HO-s,  whereas  cell  death  and  detachment  (black  areas  represent  cell-free  regions)  was  noted  for  HCFs  cultured  with  Phutopia.  For  all  images,  scale 
bat  -  200  pm  as  seen  in  boft<nn  right  Image.  (B)  Quantification  of  live  cells  elucidated  that  Photopia  counteracts  cell  proliferation  as  live  cell  count 
did  not  sign  lit  canity  change  from  day  l  to  in.  Control,  blue:  Photopia,  ted:  artificial  iris,  green  bars. 
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Fkhjm;  6,  (A)  HC  NMH  spectra  of  Photopia  and  artificial  iris  samples  immersed  in  PBS,  The  visible  peaks  are  from  mobile  fractions  of  polymers  in 

aqueous  solution.  <B)  l*C  NMK  spectra  of  supernatant  PBS  exposed  to  Photopia  and  artificial  iris  for  I  month.  Tile  J  24-ppm  peak  seen  in  Photopia 
supernatant  is  effectively  blocked  and  not  visible  in  the  artificial  iris  supernatant  sample.  <C>  Chemical  structure  of  polyethylene,  naphthopyran,  and 
PDM.S.  The  reversible  and  photoactive  C-O  bond  in  naphthopyran  is  indicated  by  a  dashed  fine. 


our  artificial  iris  approaches  such  a  design  goal  (Fig,  3A. 
artificial  iris).  Furthermore,  when  activated  by  UV  light,  the 
artificial  iris  blocks  an  additional  10%  of  green  to  orange  light 
(450-650  nm),  thereby  enhancing  its  light  sensitivity  when 
outdoors  or  in  the  presence  of  harmful  I  TV  radiation  (Fig.  5A, 
l TV-activated  artificial  iris). 

As  the  naphthopyran  side  groups  Largely  determine  its  light- 
iransniission  properties,  Photopia 's  unique  chemical  structure 
allows  for  the  increased  attenuation  of  visible  light  in  the  active 
state,22  In  addition  to  diminishing  the  amount  of  light  entering 
the  eye  based  on  wavelength,  our  artificial  iris  has  graded  light 
transmission  depending  on  the  intensity  of  the  incident  light. 
Applying  neutral -density  filters  with  optical  density  of  0,  03, 
and  0.6  to  alter  intensity  across  the  light  spectrum  showed  Lhai 
as  the  ambient  fight  brightens,  a  larger  percentage  of  that  light 
is  absorbed  (Fig.  3B).  Thus,  similar  to  the  natural  iris,  our 
implantable  artificial  iris  has  the  sensitivity  to  inhibit  high- 
intensity  fight  while  allowing  enough  fight  to  enter  the  eye  in 
dim  settings  to  maintain  vision. 

Although  the  artificial  iris  retains  visual  acuity  through  its 
clear  pupillary  region,  it  was  hypothesized  that  the  image 
through  the  photochromie  region  may  have  decreased  clarity 
due  to  the  molecular  structure  of  Photopia,  The  classic  optics 
model  of  an  objccKartificiaJ  iriS>lcnsrimage  System  was  used  to 
characterize  the  change  in  visual  acuity  through  the  peripheral 
artificial  iris.  We  found  an  approximately  25%  decrease  in  the 
baseline  intensity  (Fig.  4),  which  corroborates  the  results  of 
previous  optical  testing  and  confirmed  the  artificial  iris’s  ability 
to  attenuate  incident  bright  light.  We  also  noted  a  quantifiable 
degree  of  light  scattering  (Fig.  4D)  as  indicated  by  the 
increased  FWllM  estimates.  However,  further  testing  indicated 
that  adjusting  the  focal  length  by  150  pm  at  which  the  image 


through  the  peripheral  artificial  iris  was  taken  could  counter 
and  minimize  this  light  scattering  (Fig.  4Q  while  providing  Lhe 
expected  fight  attenuation.  The  graph  of  pixel  intensity  verses 
position  along  a  preselected  fine  (Figs.  4A-C)  more  clearly 
demonstrates  that  changing  the  local  length  can  restore  image 
resolution  (Fig.  4D).  For  this  graph,  the  image  resolution  is 
depicted  by  the  change  in  pixel  intensity  between  the  white 
background  and  the  black  lines  as  seen  in  the  control  image.  In 
comparison,  the  artificial  iris  curve  has  decreased  baseline 
pixel  intensity  and  a  more  gradual  change  in  the  slope 
representing  the  fight  diffraction  of  the  lines  imaged  through 
the  peripheral  artificial  iris.  However,  the  similarity  in 
steepness  of  the  pixel  intensity  curves  for  the  control  and 


Tame.  ‘*C  NMR  Chemical  Shift  and  Assignments 


Peak  Number 

Chemical  Shift,  ppm 

Assignments 

Artificial  iris 

V 

1.02 

CH*  (PDMS) 

1 

12,13 

CH*  (PE) 

2 

20.9 

CH2  (PE) 

3 

29.5 

CH2  (PE) 

4 

36 

CH2  (PE) 

5 

57.3 

CH-O  (NTHP) 

Photopia 

1 

14.29 

CH*  (PE) 

2 

23.10 

CMZ  (PE) 

3 

31,2fi 

CH>  CPE) 

4 

37.86 

5 

6032 

CH-O  (N TUP) 

Artificial  iris 

t 

5304 

CH-O  (NTHP) 

Photopia 

v 

5304 

CH-O  (NTHP) 

r 

124 

C  (NTHP) 
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artificial  iris  refocused  images  signifies  that  the  image 
degradation  can  be  reversed  and  corrected  with  a  change  in 
focal  length 

Collectively,  the  visual  acuity  data  Indicated  that  images 
taken  through  the  periphery  of  the  artificial  iris  would  have 
reduced  pixel  intensity  hut  diminished  resolution  compared 
with  images  through  the  clear  PDMS  center  unless  the  focal 
length  was  adjusted.  Fortunately,  simulating  a  change  in  focal 
length  (l.e.,  optical  power)  to  maintain  focus  on  objects  as 
their  distance  varies  is  an  inherent  process  of  the  human  lens, 
termed  accommodation.29  Given  that  the  optical  power  of  the 
human  lens  is  measured  in  diopters,  calculated  as  the  inverse 
of  focal  length,  the  150-j.im  change  in  focal  length  required  for 
visual  clarity  is  equal  to  a  change  in  optical  power  by  0.55 
diopters.  As  the  optical  power  of  the  human  lens  can  vary  up 
to  10  diopters  or  more  until  age  *10,  it  appears  the  natural  lens 
could  accommodate  the  light  diffraction  through  the  periph¬ 
eral  artificial  iris  postimphnt.  Thus,  individuals  with  our 
artificial  iris  could  maintain  visual  acuity  when  light  passes 
through  both  the  clear  center  and  the  light -responsive 
peripheral  region  of  the  implant.10 

Once  it  was  determined  that  our  artificial  iris  met  the 
optical  design  criteria,  we  pursued  cell  culture  experiments  to 
ascertain  the  biocompatibility  of  the  artificial  iris  in  vitro. 
Human  cornea!  fibroblasts  were  indirectly  exposed  to  samples 
of  the  artificial  iris  and  Photopia  alone  for  up  to  10  days.  Cell 
viability  was  established  when  normal  morphology  and 
proliferation  rates  were  observed  for  HCFs  exposed  io  the 
artificial  iris  compared  with  control  samples  (Fig.  5 A).  The 
biocompatibility  of  the  artificial  iris  was  then  validated  when 
subsequent  analysis  indicated  no  significant  change  in  live  cell 
count  (Fig.  5b)  for  artificial  iris  and  control  samples  at  all  study 
time  points,  In  contrast,  the  decreased  number  of  HCFs  in  the 
Phoiopia  samples  highlighted  the  potential  adverse  effects  of 
exposure  to  Phoiopia  alone.  As  the  artifieial  iris  consists  of 
Phoiopia  encased  in  PDMS,  a  polymer  whose  biocompatibility 
has  been  well  established 31  12  the  in  vitro  data  suggested  that 
the  design  of  the  artificial  iris  inhibited  cell  death  despite  the 
possible  toxicity  of  Photopia.  We  postulate  that  interaction  of 
ionic  media  with  reactive  surface  naphthopyran  molecules  in 
Photopia  might  have  caused  the  cell  death  observed  with  the 
Photopia  samples. 

based  on  the  in  vitro  data,  we  endeavored  to  determine  the 
molecular  struct  Lire  of  Phoiopia  and  any  potentially  reactive 
side  chains  that  could  cause  toxicity  in  solution.  In  addition, 
we  aimed  to  ascertain  the  stability  of  the  artificial  iris  implant 
over  time.  Commonly  used  to  determine  the  structure  and 
morphology  of  inorganic  materials,  NMR  spectroscopy  detects 
the  chemical  environment  of  atomic  nuclei  based  on  the 
differences  in  electromagnetic  radiation  they  emit  in  the 
presence  of  a  high  magnetic  field.33  The  mobile  molecular 
groups  in  both  samples  were  identified  using  l3C  NMR  scans  of 
Phoiopia  and  the  artificial  iris  in  PhS  and  are  shown  in  the 
Table.  Common  peaks  between  the  scans  (Fig.  6A)  were  used 
to  identify  the  molecular  components  of  Photopia.  The  upheld 
shift  of  the  Phoiopia- related  peaks  in  the  artificial  iris  spectra 
and  the  presence  of  the  strong  Cl  1  ^  peak  associated  with  the 
PDMS  coating  indicated  that  PDMS  served  as  an  effective 
barrier  in  preventing  the  interaction  of  Photopia  with  the  ions 
in  aqueous  solution.  Another  set  of  NMR  experiments  (Fig.  6b) 
on  PBS  exposed  to  the  artificial  iris  or  Photopia  for  I  month 
was  conducted  to  ascertain  if  molecular  leaching  would  occur 
over  time.  Nuclear  magnetic  resonance  testing  of  the 
supernatant  liquid  alone  showed  a  peak  at  124  ppm  in  the 
Phoiopia  sample  not  seen  in  the  artificial  iris  scan  that  could 
have  arisen  from  the  aromatic  carbon  group  in  the  naphihtJ- 
pyran  (Fig.  6C)  This  additional  molecular  species  in  solution 
may  explain  the  cell  death  seen  in  images  of  Photopia  at  day  10 


(see  Fig.  5 A).  Furthermore,  it  confirms  that  the  PDMS  coating 
in  the  artificial  iris  effectively  prevents  leaching  of  aromatic 
carbon  group  while  ensuring  in  vitro  cell  survival  and  stability 
of  the  implant  over  time. 

Conclusions 

Our  light-activated  artificial  iris  is  a  biocompatible  implant  thaL 
partially  mimics  the  functionality  of  the  natural  iris  while 
restoring  iris  aesthetics.  Incorporation  of  a  photo-responsive 
material  encased  in  a  PDMS  matrix  to  engineer  our  device 
provides  for  a  new  artificial  iris  design.  'Hie  dynamic  light 
attenuation  based  on  the  ambient  light  intensity  and  wave¬ 
length  GW  or  visible  light)  characteristic  of  our  artificial  iris 
improves  on  the  properties  of  commercial  artificial  iris 
implants.  The  ability  to  attenuate  up  to  60%  of  incident  light 
provides  our  artificial  iris  with  the  means  to  alleviate  or  at  least 
significantly  reduce  the  symptoms  of  iris  damage,  such  as 
photophobia,  glare,  and  haloes.  By  maintaining  visual  acuity  for 
light  traversing  both  the  central  transparent  PDMS  as  well  as 
the  light-responsive  peripheral  region,  our  implant  proves 
comparable  to  commonly  used  intraocular  lens  implants. 
Unlike  photochromie  intraocular  lenses  that  filter  out  UV  and 
blue  lights  by  absorbing  them, 3 ’  ^  the  artificial  iris  we 
developed  filters  and  uses  UV  lights  to  activate  the  photo- 
response  materials  to  dynamically  modulate  the  incident  light 
intensity.  In  addition  to  providing  enhanced  optical  properties, 
it  demonstrated  biocompatibility  in  vitro  despite  potential 
toxicity  from  Photopia  alone.  Follow-up  testing  using  13C  NMR 
served  to  identify  structural  characteristics  of  Photopia 
including  ones  that  may  lead  to  toxicity  in  solution  as  well  as 
established  that  the  PDMS  coating  in  the  artificial  iris  prevents 
leaching.  Future  work  on  artificial  iris  implantation  in  the 
rabbit  eye  would  be  necessary  to  confirm  these  findings  in 
vivo.  Although  the  final  design  of  the  artificial  iris  will  be 
optimized  based  on  preclinical  animal  model  studies,  we  have 
successfully  engineered  a  photo-responsive  and  hioconipatiblc 
artificial  iris  implant  with  dynamic  light-attenuation  properties 
that  may  provide  a  new'  treatment  option  for  patients  with 
congenital  or  trauma-induced  iris  defects. 
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